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Abstract. The paper gives a review of a very recent results related to the 
Poncelet Theorem, on the occasion of its bicentennial. We are telling the story 
of one of the most beautiful theorems of Geometry, recalling for the general 
mathematical audience the dramatic historic circumstances which led to its 
discovery, a glimpse of its intrinsic appeal, and importance of its relationship 
to the dynamics of billiards within confocal conies. We focus on the three 
main issues: A) The case of Pseudo-Euclidean spaces, presenting a recent 
notion of relativistic quadrics, and applying it to the description of periodic 
trajectories of billiards within quadrics. B) The relationship between so-called 
billiard algebra and foundations of modern discrete differential geometry which 
leads to the Double-reflection nets. C) We introduce a new class of dynamical 
systems - pseudo-integrable billiards generated by the boundary composed of 
several arcs of confocal conies having nonconvex angles. The dynamics of such 
billiards has several extraordinary properties. They are related to the interval 
exchange transformations and generate families of flows which are minimal 
but not uniquely ergodic. This type of dynamics provides a novel type of the 
Poncelet porisms — the local ones. 
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1. Introduction and history 

On November 18th, 1812, near Smolensk, a young French officer serving as a 
battery commander was wounded and his horse killed under him. This was the last 
day of the battle of Krasnoi, when the French III Corps of Marshal Ney clashed 
with three corps commanded by general Count Miloradovich of the Russian Army. 
Napoleon's army was heavily defeated and many thousands of his men were left 
during the withdrawal and imprisoned. A day after the battle, the officer was 
found by Russian soldiers. 

During his subsequent imprisonment in Saratov, which lasted from April 1813 
until June 1814, while recovering from an illness, he recalled the fundamental prin- 
ciples of geometry to which he had been introduced during his studies at Ecole 
Polytechnique. Without literature, he not only recollected what he had learned 
from his professors Monge, Carnot, and Brianchon, but also went on developing 
projective geometry, in particular properties of conies. The young officer's name 
was Jean Victor Poncelet. The notes he made in prison, called Cahiers de Sara- 
contained one of deepest, most beautiful, and most important theorems of 
projective geometry - the Great Poncelet Theorem, see [1869 . 

The first published proof of the Great Poncelet Theorem appeared in Pon- 
celet's famous work Traite des proprietes projectives des figures in 1822 [Po nl822] , 
The complete Cahiers de Saratov were published many decades later, in 1862 as 
|Ponl862"] . Meanwhile, Poncelet became a professor of Mechanics at Sorbonne and 
at the College de France, general of a brigade, the governor of the Ecole Poly- 
technique, the commander of the National Guard of the Department of the Seine, 
elected member of the Constitutional Assembly, President of the Scientific Com- 
mission for the English exhibition of 1851. Poncelet was also a Grand Officer of the 
Legion of Honour, Chevalier of the Prussian order, corresponding member of the 
academies of Sankt Petersburg, Turin, Berlin, and a Foreign member of the Royal 
Society of London. 

It is interesting to mention that, in Chebyshev's report on his business trip 
to France [Tchl852 , Poncelet is described "as a well-known scientist in practical 
mechanics" . 

Let us present now one of the formulations of the Poncelet Theorem. Suppose 
that two ellipses are given in the plane, together with a closed polygonal line in- 
scribed in one of them and circumscribed about the other one. Then, the Poncelet 



1 Saratov notebooks 
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Theorem states that infinitely many such closed polygonal lines exist — every point 
of the first ellipse is a vertex of such a polygon. Besides, all these polygons have 
the same number of sides. 

Theorem 1.1 (The Poncelet Theorem). Let C and T> be two conies in the plane. 
Suppose that there is a polygon inscribed in C and circumscribed about T>. Then 
there is infinitely many such polygons and all of them have the same number of 
sides. Moreover, each point of C is a vertex of such a polygonal line. 



Later, using the addition theorem for elliptic functions, Jacobi gave another proof 
of the theorem in 1828 (see Jacl884j). Essentially, the Poncelet Theorem is equiv- 
alent to the addition theorems for elliptic curves and Poncelet's proof represents a 
synthetic way of deriving the group structure on an elliptic curve. 

An important question is to find an analytical condition determining, for two 
given conies, if an n-polygon inscribed in one conic and circumscribed about the 
second exists. In 1853, such a condition was derived by Cayley, who used the theory 
of Abelian integrals Cayl854]. He was dealing with the Poncelet Porism in a num- 



ber of other papers [Cayl853 Ca yT855l[CayT8571|Cayl858[|Cayl861| . Cayley's work 



served as an inspiration for another great mathematician, Lebesgue, who translated 
Cayley's proof to the geometric language. He derived the proof of Cayley's condi- 
tion using methods of projective geometry and algebra, see the remarkable book 
Les coniques Lebl942]. In modern settings, Griffiths and Harris derived Cayley 
theorem by finding an analytical condition for points of finite order on an elliptic 
curve |GH1978| . 

We have to emphasize that Poncelet originally proved a statement that is much 
more general than the theorem formulated above, see [Bcrl987, Ponl822|. He de- 
rived the latter as a corollary. Namely, he considered n + 1 conies of a pencil in 
the projective plane. If there exists an n-polygon with vertices lying on the first of 
these conies and each side touching one of the other n conies, then infinitely many 
such polygons exist. We shall refer to this statement as the Full Poncelet Theorem 
and call such polygons the Poncelet polygons. 

A nice historical overview of the Poncelet Theorem, together with modern proofs 
and remarks is given in BKOR1987 . Various classical theorems of Poncelet type 
with short modern proofs are reviewed in [BB1996] . while the algebro-geometrical 
approach to families of the Poncelet polygons via modular curves is given in BM1993, 
IJakl993| . There are also two recent books on the subject |Fla2009lHTR2011j . 

The Poncelet Theorem has an important mechanical interpretation. The ellip- 
tical billiard [KT1991, Koz2003 ( Tab2005 is a dynamical system where a material 
point of the unit mass is moving under inertia, or in other words, with a constant 
velocity inside an ellipse and obeying the reflection law at the boundary, i.e. hav- 
ing congruent impact and reflection angles with the tangent line to the ellipse at 
any bouncing point. It is also assumed that the reflection is absolutely elastic and 
friction is neglected. 

It is well known that any segment of a given elliptical billiard trajectory is tangent 
to the same conic, confocal with the boundary CCS1 993] . If a trajectory becomes 
closed after n reflections, then the Poncelet Theorem implies that any trajectory of 
the billiard system, which shares the same caustic curve, is also periodic with the 
period n. 
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The Full Poncelet Theorem also has a mechanical meaning. The configuration 
dual to a pencil of conies in the plane is a family of confocal second order curves 
Arnl978 . Let us consider the following, slightly unusual billiard. Suppose n con- 
focal conies are given. A particle is bouncing on each of these n conies respectively. 
Any segment of such a trajectory is tangent to the same conic confocal with the 
given n curves. If the motion becomes closed after n reflections, then, by the Full 
Poncelet Theorem, any such a trajectory with the same caustic is also closed. 

The statement dual to the Full Poncelet Theorem can be generalized to the 
(/-dimensional space [CCS1993] (see also |Prel999|IPrej ). Suppose vertices of the 
polygon X1X2 ■ ■ .x n are respectively placed on confocal quadric hyper-surfaces Qi, 
Q2, • • • j Qn in the d-dimensional Euclidean space, with consecutive sides obeying 
the reflection law at the corresponding hyper-surface. Then all sides are tangent to 
some quadrics Q 1 , . . . , confocal with {Qi}; for the hyper-surfaces {Qi, Q J '}, 

an infinite family of polygons with the same properties exist. Systematic exposition 
of this higher-dimensional theory has been presented in a survey paper [DR2010 
and in the book [DR2011) . 

But, more than one century before these quite recent results, in 1870, Darboux 
proved the generalization of the Poncelet Theorem for a billiard within an ellipsoid 
in the three-dimensional space Darl870 . It seems that his work on this topic was 
completely forgotten until very recently, see |DR201l] . 

It is natural to search for a Cayley-type condition related to generalizations of the 
Poncelet Theorem. Such conditions for the billiard system inside an ellipsoid in the 
Eucledean space of arbitrary finite dimension were derive d in [DR1998a ( DR1998b]. 
In recent papers lDR2004ilDR2005llDR20Q6bllDR2006aj . algebro-geometric condi- 
tions for existence of periodical billiard trajectories within k quadrics in (/-dimen- 
sional Euclidean space were derived. 

Most of the results on the subject obtained by 2008 have been presented in the 
book |DR201l] . An important part of the book has been devoted to the Griffiths- 
Harris program of development of a synthetic approach to higher genera addition 
theorems. The book also offered a historical overview of the subject, and it included 
a detailed analysis of the results of Darboux and the contributions of Lebesgue. 

About the same time, an extremely interesting book |l)ui20H)i appeared, devoted 
to discrete integrable systems from the point of view of the QRT mapto, elliptic 
surfaces, and elliptic fibrations. That book provides approach to the Poncelet 
Theorem as to an important example of the symmetric QRT maps, see |Dui2010] 
Chapter 10]. Connections to the elliptic billiards are presented in [Dui2010 l Chapter 
11.2]. 

The present paper is devoted to the bicentennial jubilee of the celebrated Pon- 
celet Theorem and it mostly exposes the current advances of the subject, the results 
which have been obtained in the last four years. The interrelation between geom- 
etry of pencils of quadrics and related billiard dynamics is continuing to play a 
crucial role. 

The projective geometry nucleus of that billiard dynamics is the Double reflection 
theorem, see Theorem 12.91 below. There are four lines belonging to a certain linear 
space and forming a double reflection configuration: these four lines reflect to each 
other according to the billiard law at some confocal quadrics, see Section l2~2l 



2 The QRT maps are named after Quispel, Roberts, and Thompson QRT1988 . 
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The double reflection configuration is a cornerstone of a new type of integrable 
line congruences - so-called double reflection nets. In these discrete integrable 
systems of geometric origin, double reflection configurations are playing the role 
of quad-equations. The integrability condition is a consequence of an operational 
consistency of the billiard algebra from [DR2008 , see the Six-pointed star theorem 
- Theorem 14.11 These results are presented in Section |4] 

Section |3] is devoted to pencils of quadrics in the pseudo-Euclidean spaces and to 
related billiard dynamics and Poncelet configurations. A novelty of our approach 
is based on the notion of relativistic quadrics, see [DR2012a] , The general the- 
ory of relativistic quadrics is exposed in the Section 13.21 Section 13.11 is devoted 
to the case of the Minkowski plane and serves as an introduction to the higher- 
dimensional cases. It also contains a detailed and quite elementary description of 
periodic trajectories of elliptic billiards and an analysis of topological properties of 
such systems. Generalized Cayley-type conditions for pseudo-Euclidean spaces of 
arbitrary dimension are derived in Section 13.31 

The last section is devoted to the billiards in Euclidean plane with more complex 
boundary, formed by arcs of conies from a confocal family and related Poncelet 
Porism, see [DR2012cj . 

For the case of billiard systems within confocal conies without non-convex angles 
on the boundary, it is well known that the famous Poncelet Porism holds: 

(A) if there is a periodic billiard trajectory with one initial point of the boun- 
dary, then there are infinitely many such periodic trajectories with the same 
period, sharing the same caustic; 

(B) even more is true, if there is one periodic trajectory, then all trajectories 
sharing the same caustic are periodic with the same period. 

See |DR2006al[DR2006b[|DRlOTT] where the corresponding conditions of Cayley- 
type were derived. 

However, when non-convex angles exist on the boundary, which is the case stud- 
ied in Section O one sees that (A) is still generally true. However, (B) is not true 
any more. The Poncelet Porism is true locally, but not globally. Algebro-geometric 
conditions of Cayley's type in such a case provide only sufficient but not necessary 
conditions for periodicity. A deeper analysis of dynamics in this case is related to 
the class of interval exchange transformations and to the use of a modified Keane's 
condition. Section 15.61 is concluded by an explicit example of a billiard system 
which satisfies the Cayley type condition, but is still not periodic since it satisfies 
the Keane's condition as well. In Section r5.71 we show that there are infinitely many 
billiard tables bounded by arcs of confocal conies such that the corresponding flow 
will not be uniquely ergodic. In particular, Theorem 15.161 shows that the interval 
exchange transformation corresponding to certain tables is equivalent to the Veech 
example of minimal and not uniquely ergodic systems [Veel969 ( MT2002 . 

As a starting point, we collect some basic notions and facts, which are going to 
be used in sequel. 



2. Billiards and quadrics 
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2.1. Elliptical billiards and confocal conies in the Euclidean plane. We 

will start with the billiard systems within confocal conies as a mechanical setting 
for the Poncelet porism. 

A general family of confocal conies in the plane can be represented in the fol- 
lowing way: 

(2.1) C x : + JlL. = 1, A e R, 

a — A o — A 

with a > b > being constants. 

By the famous Chasles' theorem Chal827], each segment of a given billiard 
trajectory is tangent to a fixed conic that is confocal to the boundary (see also 
[KT1 991 . DR 2011j ). This conic a caustic of the given trajectory. 

Now, fix a constant ao < b and consider billiard trajectories within confocal 
ellipses C\ (A < ao) having ellipse C ao as caustic. 

Proposition 2.1. There exist metrics p on a conic C ao and a function 

p : (— oo, «o) R 

satisfying: 

• M(C ao ) = !; 

• metric p is non-atomic, i.e. p({X}) = for each point X on C ao ; 

• n(t) for each open arc I of C aa ; 

• for any A < a , and each triplet of points X g C ao , Y G C ao , A £ C\. such 
that segments XA and AY satisfy the reflection law on C\, the following 
equality holds: 

H(XY) = p(X). 

Proof. Take Ao such that there is a closed billiard trajectory in C\ with caustic 
Ca ■ By |Kinl994j , there is a metric fx satisfying the requested properties for A = Ao 
- moreover, such a metric is unique up to multiplication by a constant. By this 
uniqueness property and Darboux theorem on grids |Darl914] (see also DR2006b ( 
DR2008, DR2011J), it follows that metric fi satisfies the properties for each C\ 
having closed billiard trajectories with caustic C ao . 

For a periodic trajectory which becomes closed after n bounces on C\ and m 

vn 

windings about C\ , p(X) = — . Since rational numbers are dense in the reals, p, 

n 

will have the required properties for all A < a - D 

Remark 2.2. The function p from Provosition \2. 1\ is called the rotation function 
and its values the rotation numbers. Note that p is a continuously strictly decreasing 

function with ^0, — ^ as image: 

lim p(X) = -, lim p(A) = 0. 

A— V — oo Z A— VQq 

Elliptical billiard as a Hamiltonian system. The standard Poisson bracket for the 
billiard system is defined as: 

If g } = ^L^l _ ^£^1 + ^f_^9_ _ df_dg_ 

' dx dx dx dx dy dy dy dy 
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Define the following functions: 

x 2 y 2 (xy — yxY 



K x (x,y,x,y) 



a - A b- A (a-A)(&-A)' 

These functions represent well-known first integrals of the billiard systems, see 
[KT1991] . 

Proposition 2.3. Each two functions K\ commute: 

{K Xl ,K X2 } = 
and for Ai =/= A 2 , they are functionally independent. 

It is straightforward to prove the following 

Proposition 2.4. Along a billiard trajectory within any conic C\ , with caustic 
C ao and the speed of the billiard particle being equal to s, the value of each function 
K\ is constant and equal to 

K = U ° ~ A 2 
A (a- X)(b- A) ' ' 

Corollary 2.5. Each K\ is integral for the billiard motion in any domain with 
border composed of a few arcs of confocal conies. 

2.2. Double reflection configurations. Billiards within pencils of quadrics in- 
duce fruitful dynamical systems in arbitrary dimension. They are meaningful in 
spaces with non-Euclidean metric as well, and even in spaces without any metric 
at all. 

In this section, we review a fundamental projective geometry configuration of 
double reflection in the d-dimensional projective space P d over an arbitrary field of 
characteristic not equal to 2. A detailed discussion on this matter can be found in 
|DR2008l[DR20TT] (see also |CCS1993] 1. 

The section is concluded by Proposition ^. 141 where we show that double reflec- 
tion configuration can take the role of the quad-equation, that is every line in such 
a configuration is determined by the remaining three. This simple observation is 
going to play a key role in Section 01 in particular in Section [4.21 

Let us start with recalling the notions of quadrics and confocal families in the 
projective space. 

A quadric in P'' is the set given by equation of the form: 

where Q is a symmetric (<f + 1) X (d + 1) matrix, and £ = [£o '■ £1 '■ • • • '■ £d] are 
homogeneous coordinates of a point in the space. 
Assume two quadrics are given: 

Si : (Qi£,£)=0, Q 2 : (Q 2 £,0 = 0. 

A pencil of quadrics is the family of quadrics given by equations: 

((Qi + AQ a )£,O = 0, AeP 1 . 

A confocal system of quadrics is a family of quadrics such that its projective dual 
is a pencil of quadrics. 
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Now, let us recall the definition of reflection off a quadric in the projective space, 
where metrics is not defined. This definition, together with its crucial properties 
- the One reflection theorem and the Double reflection theorem, can be found in 
[UUS1993] . 

Denote by u the tangent plane to Q\ at point x and by z the pole of u with 
respect to Q.2- Suppose lines £\ and £2 intersect at x, and the plane containing 
these two lines meet u along I. 

Definition 2.6. If lines £\, £2, xz, £ are coplanar and harmonically conjugated, 
we say that £\ is reflected to £2 off quadric Q\. 

It can be proved that this definition does not depend on the choice of quadric 
Q2 from a given confocal system [CCS1993] . 

If we introduce a coordinate system in which quadrics Q\ and Q2 are confocal 
in the usual sense, the reflection introduced by Definition 12.61 is the same as the 
standard, metric one. 

Theorem 2.7 (One reflection theorem). Suppose line £\ is reflected to £2 off Q\ 
at point x, with respect to the confocal system determined by quadrics Q\ and Q.2- 
Let £\ intersect Q2 at y[ and y±, u be the tangent plane to Q\ at x, and z the pole 
of u with respect to Q.2- Then lines y[z and y\z respectively contain intersecting 
points y' 2 and 7/2 of line £2 with Q2 . The converse is also true. 

Corollary 2.8. Let lines £\ and £2 reflect to each other off Qi with respect to the 
confocal system determined by quadrics Qi and Q.2- Then £\ is tangent to Q2 if 
and only if £2 is tangent to Q2; £\ intersects Q2 at two points if and only if £2 
intersects Q2 at two points. 

Next theorem is crucial for our further considerations. 

Theorem 2.9 (Double reflection theorem) . Suppose that Qi, Q2 are given quadrics 
and X\Q\, yi € Q.2- Let u\ be the tangent plane of Qi at x\; z\ the pole of u\ with 
respect to Q2; V\ the tangent plane of Q2 at y\; and w± the pole of Vi with respect 
to Qx- Denote by X2 the intersecting point of line W1X1 with Q\, X2 7^ Xi; by y2 
the intersection of y\Z\ with Q2, y2 7^ yi; and £\ — x%yi, £2 — X\y2, £\ = y\X2, 

£'2 = X2V2- 

Then pair £\, £2 obey the reflection law off Q\ at x\; £\, i\ obey the reflection 
law off Q2 at y\; £2, £' 2 obey the reflection law off Q2 at 2/2/ and £\, £' 2 obey the 
reflection law off Qi at point X2 ■ 

Let us remark that in Theorem 12.91 the four tangent planes at the reflection 
points belong to a pencil. 

Corollary 2.10. // line £\ is tangent to a quadric Q! confocal with Qi and Q2, 
then £2, £'\, £'2 also touch Q! . 

The following definition of virtual reflection configuration and double reflection 
configuration is from [DR2008 . These configurations played the central role in 
that work. In Thcorcm l2.121 some important properties of these configurations are 
given. This theorem is proved in [DR2008] . 

The configurations are connected with so-called real and virtual reflections, but 
their properties remain in the projective case, too. 

Let points x\, X2 belong to Qi and yi, ?/2 to Q2. 
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Definition 2.11. We will say that the quadruple of points Xi,X2,yi,y2 constitutes 
a virtual reflection configuration if pairs of lines x\y\, x\y<i; x%yx, xiyi; x\y\, 
Xiyi' } x\y%, x%y% satisfy the reflection law at points x±, X2 off Q\ and y\, y% off Q2 
respectively, with respect to the confocal system determined by Q\ and Q.2- 

If, additionally, the tangent planes to Q±, Q2 at x\, Xi; 2/1,2/2 belong to a pencil, 
we say that these points constitute a double reflection configuration (see Figure^). 




Figure 1. Double reflection configuration 

Now, we list some of the basic facts about double reflection configurations. 

Theorem 2.12. Let Qi, Q2 be two quadrics in the projective space P d , x\, X2 
points on Q\ and y\, 2/2 on Q.2- If the tangent hyperplanes at these points to 
the quadrics belong to a pencil, then x\,X2,yi,y2 constitute a virtual reflection 
configuration. 

Furthermore, suppose that the projective space is defined over the field of reals. 
Introduce a coordinate system, such that Q\, Q2 become confocal ellipsoids in the 
Euclidean space. If Q\ is placed inside Q2, then the sides of the quadrilateral 
X\y\X2V2 obey the real reflection from Q2 and the virtual reflection from Q\. 

The statement converse to Theorem 12.121 is the following 

Proposition 2.13. In the Euclidean space ~E d , two confocal ellipsoids 8\ and 82 
are given. Let points X\, X2 belong to £%, Y\, Y% to £2, and let ai, a.i, f3i, P2 be the 
corresponding tangent planes. If a quadruple Xi, X2, Y"i, I2 is a virtual reflection 
configuration, then planes ct\, 0:2, Pi, /?2 belong to a pencil. 

The next proposition shows that three lines of a double reflection configuration 
uniquely determine the fourth one. 

Proposition 2.14. Let I, i\, £2 be lines and Q\, Q2 quadrics in the projective 
space. Suppose that I, l\ reflect to each other off Qi, and £, £2 off Q2, with respect 
to the confocal system determined by these two quadrics. Then there is a unique 
line £12 such that four lines I, l\, £2, £12 form a double reflection configuration. 

Remark 2.15. Proposition \2.ll\ shows that double reflection configuration is play- 
ing the role of the quad-equation for lines in the projective space. 
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2.3. Pseudo-Euclidean spaces and confocal families of quadrics. In this 
section, we first give a necessary account of basic notions connected with pseudo- 
Euclidean spaces and their confocal families of quadrics. 

Pseudo-Euclidean spaces. The pseudo- Euclidean space E fc,z is a d-dimensional space 
R d with the pseudo- Euclidean scalar product: 

(2.2) (a;, y) k ,i = X\V\ H V x k y k - x k+1 y k+ i x d y d . 

Here, k, I G {1, . . . , d — 1}, k + I — d. The pair (fc, I) is called the signature of the 
space. Denote E kd = diag(l, 1, . . . , 1, — 1, . . . , —1), with k l's and I —Vs. Then the 
pseudo-Euclidean scalar product is: 

(x,y)k,l = E kA x o y, 

where o is the standard Euclidean product. 

The pseudo-Euclidean distance between points x, y is: 



Since the scalar product can be negative, notice that the pseudo-Euclidean distance 
can have imaginary values as well. 

Let £ be a line in the pseudo-Euclidean space, and v its vector. £ is called: 

• space-like if (v,v) k .i > 0; 

• time-like if (v,v)k,i < 0; 

• and light-like if (v,v) k j — 0. 

Two vectors x, y are orthogonal in the pseudo-Euclidean space if (x, y)k,l — 0. Note 
that a light-like line is orthogonal to itself. 

For a given vector » / 0, consider a hyper-plane vox = 0. Vector E kd v is 
orthogonal to the hyper-plane; moreover, all other orthogonal vectors are collinear 
with Ek,iv. If v is light-like, then so is E k .iv, and Ekjv belongs to the hyper-plane. 

Billiard reflection in the pseudo -Euclidean space. Let v be a vector and a a hyper- 
plane in the pseudo-Euclidean space. Decompose vector v into the sum v = a + n a 
of a vector n a orthogonal to a and a belonging to a. Then vector v' = a — n a is 
the billiard reflection of n on a. It is easy to see that then v is also the billiard 
reflection of v' with respect to a. 

Moreover, let us note that lines containing vectores v, v' , a, n a are harmonically 
conjugated 1KT2009] . 

Note that v = v' if v is contained in a and v' = —v if it is orthogonal to a. If n a 
is light-like, which means that it belongs to a, then the reflection is not defined. 

Line £' is the billiard reflection of £ off a smooth surface S if their intersection 
point £ n £' belongs to S and the vectors of £, £' are reflections of each other with 
respect to the tangent plane of S at this point. 

Remark 2.16. It can be seen directly from the definition of reflection that the type 
of line is preserved by the billiard reflection. Thus, the lines containing segments 
of a given billiard trajectory within S are all of the same type: they are all either 
space-like, time-like, or light-like. 
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If S is an ellipsoid, then it is possible to extend the reflection mapping to those 
points where the tangent planes contain the orthogonal vectors. At such points, a 
vector reflects into the opposite one, i.e. v' — —v and £' — £. For the explanation, 
see |KT2009] . As follows from the explanation given there, it is natural to consider 
each such reflection as two reflections. 

Families of confocal quadrics. For a given set of positive constants 01, a 2 , . . . , 
ad, an ellipsoid is given by: 

22 2 

(2.3) £ : a + £2 + ... + £d =L 

oi a 2 ad 

Let us remark that equation of any ellipsoid in the pseudo-Euclidean space can be 
brought into the canonical form (12.3[) using transformations that preserve the scalar 
product (|2~2"j) . 

The family of quadrics confocal with £ is: 

(2.4) Qx : + + + ■ ■ ■ H ^L— = i ; A e R. 

ai - A afc — A afc + i + A ad + A 

Unless stated differently, we are going to consider the non-degenerate case, when 
set {ai, . . . , afc, — dfc+i, . . . , —ad} consists of d different values: 

ai > a 2 > • • • > a k > > -Ofe+i > • • • > -a d . 

For A € {ai, . . . , a^, — afc+i, . . . , — ad}, the quadric Qa is degenerate and it coin- 
cides with the corresponding coordinate hyper-plane. 

It is natural to join one more degenerate quadric to the family (|2.4[) : the one 
corresponding to the value A = 00, that is the hyper-plane at the infinity. 

For each point x in the space, there are exactly d values of A, such that the 
relation (|2.4p is satisfied. However, not all the values are necessarily real: either all 
d of them are real or there are d — 2 real and 2 conjugate complex values. Thus, 
through every point in the space, there are either d or d — 2 quadrics from the 
family (H3J) [KT2009] . 

The line x + tv (t € R) is tangent to quadric Qx if quadratic equation: 

(2.5) A x (x + tv) o (x + tv) = 1, 

has a double root. Here we denoted: 

i ( 1 11 1 
A x = diag -, • • • , -, — , • • • , — 

\ai— A ak - A a k+ i + A a d + A 
Now, calculating the discriminant of (|2.5p . we get: 

(2.6) (^ A a; o v) 2 - {Axvov){A x xox - 1) = 0, 
which is equivalent to: 

d 



(2.7) Y, \ =0, 



=1 

where 

(2.8) F &>")= ^ + Y, { 7 3 a~-ef 
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with e's given by: 




The equation (|2.7|) can be transformed to: 
(2.9) d V{X) =0, 

where the coefficient of A d_1 in V(X) is equal to (v,v)f-i- Thus, polynomial V(X) is 
of degree d — 1 for space-like and time-like lines, and of a smaller degree for light- 
like lines. However, in the latter case, it turns out to be natural to consider the 
polynomial V{\) also as of degree d—1, taking the corresponding roots to be equal 
to infinity. So, light-like lines are characterized by being tangent to the quadric 

Qoo- 

Having this setting in mind, we note that it is proved in KT2009 that the 
polynomial V(X) has at least d — 3 roots in R U {oo}. 
Thus, we have: 

Proposition 2.17. Any line in the space is tangent to either d—1 ord—3 quadrics 
of the family ^-4\ l- U this number is equal to d — 3, then there are two conjugate 
complex values of X, such that the line is tangent also to these two quadrics in C d . 

This statement with the proof is given in [KT2009] . Let us remark that in 
KT2009 is claimed that light-like line have only d — 2 or d — 4 caustic quadrics. 
That is because Qoo is not considered there as a member of the confocal family. 

As noted in [KT2009 , a line having non-empty intersection with an ellipsoid 
from (|2.4p will be tangent to d — 1 quadrics from the confocal family. However, we 
are going to prove this in another way, which will provide us some more insight into 
the distribution of the parameters of the caustics along the real axis. Next theorem 
will also contain a detailed description of distribution of parameters of quadrics 
containing given point placed inside an ellipsoid from (j2,4[) . 

Theorem 2.18. In pseudo-Euclidean space E ' consider a line intersecting ellip- 
soid £ i2.3\) . Then this line is touching d—1 quadrics from {2.4]) - If we denote 
their parameters by a,\, . . . , a d -i and take: 

{bi, bp, ci, c q } = {eiai, s d a d , ati, a d -i}, 

c q < ■ ■ ■ < c 2 < ci < < bi < b 2 < ■ ■ ■ < b p , p + q = 2d— 1, 

we will additionally have: 

• if the line is space-like, then p = 2k — 1, q = 21, ai = b p , on G {&2j-ij &2i} 
for 1 < i < k — 1, and otj+k-i € {cjy-i, C2j} for 1 < j < I; 

• if the line is time-like, then p = 2k, q = 21 — 1, c q — —a d , ai € {£>2i-i, &2i} 
for 1 < i < k, and Oj+fc € {c 2 j-i, C2j} for 1 < j < I — 1; 

• if the line is light-like, then p — 2k, q = 21 — 1, b p = oo = a^, = a\, 
on € {b 2l -i,b 2 i} fori < i < k-1, anda J+k G {c 2j -i, c 2j } for 1 <j< l-l. 

Moreover, for each point on £ inside £, there is exactly d distinct quadrics from 
\2.J$ containing it. More precisely, there is exactly one parameter of these quadrics 
in each of the intervals: 

(c 2 /_i,C2/_ 2 ), (c 3 ,c 2 ), (ci,0), (0, 6i), (b 2 ,b 3 ), (b 2 k-2, hk-i)- 
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The analogue of Theorem 12.181 for the Euclidean space, is proved in |Audl994] . 

Corollary 2.19. For each point placed inside an ellipsoid in the pseudo-Euclidean 
space, there are exactly two other ellipsoids from the confocal family containing this 
point. 

3. Pseudo-Euclidean spaces and Poncelet theorem 
3.1. Minkowski plane, confocal conies and billiards. 

Confocal conies in the Minkowski plane. Here, we give a review of basic properties 
of families of confocal conies in the Minkowski plane, sec [DR2012a . 
Denote by 

(3.1) S:^ + £ = l 

a b 

an ellipse in the plane, with a, b being fixed positive numbers. 
The associated family of confocal conies is: 

(3.2) C X : + = 1, A e R. 

a — A o + A 

The family is shown on Figure [2] We may distinguish the following three sub- 



FiGURE 2. Family of confocal conies in the Minkowski plane. 

families in the family p.2[) : 

• for A e (—6, a), conic C\ is an ellipse; 

• for A < —b, conic C\ is a hyperbola with a;-axis as the major one; 

• for A > a, it is a hyperbola again, but now its major axis is y-axis. 
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In addition, there are three degenerated quadrics: C a , Cb, Coo corresponding to 
y-axis, x-axis, and the line at the infinity respectively. Note the following three 
pairs of foci: F^y/a + b, 0), F 2 (-y/a + 6,0); Gi(0, \Ja + b), G 2 (0, -\[a + b); and 
Hi(l : -1 : 0), H 2 (l : 1 : 0) on the line at the infinity. 
We notice four distinguished lines: 

x + y = va + 6, x + y = — va + 6, 

x — y = \J a + b, x — y = -va + b. 

These lines are common tangents to all conies from the confocal family. 

Conies in the Minkowski plane have geometric properties analogous to the conies 
in the Euclidean plane. Namely, for each point on conic C\, either sum or difference 
of its Minkowski distances from the foci F\ and F 2 is equal to 2y 'a — A; either sum 
or difference of the distances from the other pair of foci G±, G 2 is equal to 2^/—b — A 
|DR2012aj . 

Relativistic conies. Since a family of confocal conies in the Minkowski plane con- 
tains three geometric types of conies, it is natural to introduce relativistic conies, 
which arc suggested in BM1962 . In this section, we give a brief account of related 
analysis. 

Consider points Fi(x/a + b, 0) and F 2 (—\/a + b, 0) in the plane. 
For a given constant c G R + U iR + , a relativistic ellipse is the set of points X 
satisfying: 

dist 1 , 1 (F 1 ,X) + dist ljl (F 2 ,X) = 2c, 
while a relativistic hyperbola is the union of the sets given by the following equations: 

disti,i(Fi,X) -disti,i(F 2 ,X) = 2c, 
disti,i(F 2 ,X) -dist M (Fi,X) = 2c. 

Relativistic conies can be described as follows. 

< c < \/ a + b: The corresponding relativistic conies lie on ellipse C a _ c 2 from 
family (|3.2[) . The ellipse C a „ c 2 is split into four arcs by touching points with 
the four common tangent lines; thus, the relativistic ellipse is the union of 
the two arcs intersecting the y-axis, while the relativistic hyperbola is the 
union of the other two arcs. 

c > ya + b: The relativistic conies lie on C a _ c 2 - a hyperbola with a;-axis as 
the major one. Each branch of the hyperbola is split into three arcs by 
touching points with the common tangents; thus, the relativistic ellipse is 
the union of the two finite arcs, while the relativistic hyperbola is the union 
of the four infinite ones. 

c is imaginary: The relativistic conies lie on hyperbola C a _ c 2 - a hyperbola 
with y-axis as the major one. As in the previous case, the branches are 
split into six arcs in total by common points with the four tangents. The 
relativistic ellipse is the union of the four infinite arcs, while the relativistic 
hyperbola is the union of the two finite ones. 

The conies are shown on Figure [3l 

Notice that all relativistic ellipses are disjoint with each other, as well as all 
relativistic hyperbolas. Moreover, at the intersection point of a relativistic ellipse 
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Figure 3. Relativistic conies in the Minkowski plane: relativistic 
ellipses are represented by full lines, and hyperbolas by dashed 
ones. 

which is a part of the geometric conic Ca x from the confocal family (|3.2[) and a 
relativistic hyperbola belonging to C\ 2 , it is always Ai < A2. 

Periodic trajectories of elliptical billiard. Analytic conditions for existence of closed 
polygonal lines inscribed in one conic and circumscribed about another one in the 
projective plane are derived by Cayley [Cayl854, Cayl86i"] . They can be applied to 
billiard trajectories within ellipses in the Minkowski plane as well, since each such 
trajectory has a caustic among confocal conies. In this section, we shall analyze in 
more detail some particular properties related to the Minkowski geometry. 

Theorem 3.1. In the Minkowski plane, consider a billiard trajectory T within 
ellipse £ given by equation US. 

The trajectory is periodic with period n = 2m if and only if the following condi- 
tion is satisfied: 







( B 3 


B 4 


B m +i ^ 




(3.3) 


det 




B 5 


B m +2 


= 






\ 


B m +2 


B2m-\ ) 





Trajectory T is periodic with period n = 2m + 1 if and only if C a is an ellipse and 
the following condition is satisfied: 







/ £3 


B 4 


B m +2 ^ 




(3.4) 


det 


B m +i 

\ Cm+l 


B m +2 
C m +2 


B2m 
C2m J 


= 
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Here: 

y/{a-t)(b + t)(a-t) = B + B x t + B 2 t 2 + ..., 



(°-W + «) =Cb + Cl t + C ^ + ... 

a — t 

are the Taylor expansions around A = 0. 

Proof. Each point inside £ is the intersection of exactly two ellipses and C\ 2 
from (|3.2j) . Parameters Ai, A2 are generalized Jacobi coordinates. Take Ai < A2. 

Consider first the case when C a is a hyperbola. Then along T these coordinates 
will take values in segments [—6, 0] and [0, a] respectively with the endpoints of the 
segments as the only local extrema. Ai achieves value —b at the intersections of 
T with the x-axis, while A2 achieves a at the intersections with y-axis. At each 
reflection point, one of the coordinates achieves value 0. They can both be equal to 
only at the points where £ has a light-like tangent, and there reflection is counted 
twice. 

This means that on a closed trajectory the number of reflections is equal to the 
number of intersection points with the coordinate axes. Since a periodic trajectory 
crosses each of the coordinate axes even number of times, the first part of the 
theorem is proved. 

The condition on T to become closed after n reflections on £ , n\ crossings over 
x-axis, and n 2 over y-axis is that the following equality 

niP a + n 2 P-b = nP a 

holds on the elliptic curve: 

s 2 = (a-t)(b + t)(a-t), 

where by Pp we denoted a point on the curve correspondig to t = /3, and Px, is 
taken to be the neutral for the elliptic curve group. 

From the previous discussion, ni+n 2 = n and all three numbers are even. P a and 
P-b are branching points of the curve, thus 2P a = 2P_/, = 2P OQl so the condition 
becomes nP = nP^, which is equivalent to (|3.3p . 

Now suppose C a is an ellipse. The generalized Jacobi coordinates take values in 
segments [—6,0], [0, a] or in [a,0], [0,a], depending on the sign of a. Since both 
cases are proceeded in a similar way, we assume a < 0. 

Coordinate Ai has extrema on T at the touching points with the caustic and 
some of the reflection points, while A2 has extrema at the crossing points with y-axis 
and some of the reflection points. 

The condition on T to become closed after n reflections on £ , with n± crossings 
over y-axis, and n 2 touching points with the caustic is: 

niP a + n 2 P a = nP Ql 

with n\ + n 2 = n and n\ even. 

Thus, for n even we get (|3 . 3|) in the same manner as for a hyperbola as a caustic. 

For n odd, the condition is equivalent to nPo — (n — \)Poo + P a - Notice that 
one basis of the space £((n— l)Poo + P a ) is: 

It t m s ts t m ~ 2 s S 

5 5''*) ) °! vo i * ' " ) L °1 1 

t — a 
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Using this basis, as it is shown in [DR2006a, DR2011]. we obtain (|3.4|) . □ 

Example 3.2 (3-periodic trajectories). Let us find all Z-periodic trajectories within 
ellipse £ given by \3.1)) in the Minkowski plane, i. e. all conies C a from the confocal 
family i3.S\) corresponding to such trajectories. 
The condition is: 

_ 3a 2 b 2 + 2a 2 ba - 2ab 2 a - a 2 a 2 - 2aba 2 - b 2 a 2 _ 
2 ~ 8{abf/ 2 a 5 / 2 ~ ' 

which gives the following solutions for the parameter a of the caustic: 

ab 



a 2 



(a + b) 2 
ab 



(a-b- 2Va 2 + ab + b 2 ), 



(a~b + 2y/a 2 +ab + b 2 ). 



(a + b) 2 

Notice that —b < a\ < < ai < a so both caustics C ai , C a2 are ellipses. 

Example 3.3 (4-periodic trajectories). By Theorem VS.lK the condition is = 0. 
Since 

(—ab — aa + ba)(—ab + aa + ba)(ab + aa + ba) 
3 = 16(a6a) 5 / 2 ' 

we obtain the following solutions: 

ab ab ab 



"l = i j «2 = — — r, a 3 = 



b — a a + b ' a + b 

Since ot\ (—6, a) and 0.2,0.3 € (—6, a), conic C ai is a hyperbola, while C a . 2 , C a3 
are ellipses. 

Example 3.4 (5-periodic trajectories). The condition is: 

?:)=»• 

Taking a = b = 1, we get that this is equivalent to 64a 6 — 16a 4 — 52a 2 + 5 = 0. This 
equation has four solutions in R, all four contained in (—1, 1), and two conjugated 
solutions in C. 

Light-like trajectories of elliptical billiard. In this section we consider in more detail 
light-like trajectories of elliptical billiard, see Figure |4] for an example of such a 
trajectory. We are going to review results from [DR2012a and illustrate them by 




Figure 4. Light-like billiard trajectory. 
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some examples. 

Periodic light-like trajectories. Let us first notice that successive segments of light- 
like billiard trajectories are orthogonal to each other (see Figure|3|, thus a trajectory 
can close only after even number of reflections. 

The analytic condition for n-periodicity of light-like billiard trajectory within 
the ellipse £ given by equation (|3.1[) can be derived as in Theorem 13.11 We get 
the condition stated in (|3.3p . with a — oo, i.e. (Bi) are coefficients in the Taylor 
expansion around t = of y/Ja — t)(b + t) = Bq + Bit + B^t 2 + . . . . 

Now, we are going to derive analytic condition for periodic light-like trajectories 
in another way, which will lead to a more compact form of i 



Theorem 3.5. Light-like billiard trajectory within ellipse £ is periodic with period 
n, where n is an even integer if and only if 

[a ( kir , n t, n\ 1 

(3.5) arctan^-ej- !<*<-, (*,-)=! j. 

As an immediate consequence, we get 

Corollary 3.6. For a given even integer n, the number of different ratios of the 
axes of ellipses having n-periodic light-like billiard trajectories is equal to: 

{ip(n)/2 if n is not divisible by 4, 
ip(n)/A if n is divisible by 4. 

ip is the Euler's totient function, i.e. the number of positive integers not exceeding 
n that are relatively prime to n. 

Remark 3.7. There are 4 points on £ where the tangents are light-like. Those 
points cut 4 arcs on £ . An n-periodic trajectory within £ hits each one of a pair of 
opposite arcs exactly k times, and — — k times the arcs from the other pair. 

Example 3.8 (4-perodic light-like trajectories). It is elementary to see that only 
circles allow such trajectory. However, we can also deduce it from condition H3.3\) . 
For n = 4, this condition reads: 

(a-b)(a + b) 2 

B 3 = 0, with B 3 ^- 



16(a&) 5 / 2 
Thus, it is equivalent to a = b. 

On the other hand, condition \3.5\) states: 



la 7r 
— = tan — = 1 . 
b 4 



Example 3.9 (6-perodic light-like trajectory). For n — 6, the condition 



which is equivalent to: 



«1H ( I" ) = I). 

i>4 -t>5 



(a - 36)(3a - b)(a + bf _ () 



{Aaby 



BICENTENNIAL OF THE GREAT PONCELET THEOREM 



19 



Thus, we get that light-like billiard trajectories are 6-periodic in ellipses with the 
ratio of the axes equal to y/3. 

From condition \S.5\l , we get the same: 



vf e { tan I' tan I} = H'7i 

A few 6-periodic trajectories are shown on Figured 




Figure 5. Light- like billiard trajectories with period 6 in the el- 
lipse satisfying a = 36. 



Example 3.10 (8-perodic light-like trajectories). For n = 8, the condition {3. 
is: 

,ui 1 1 1 )-(- t »c.t t) . u .y-j at+ " ) - ll . 



B S B e B, ) 2»(o»)»Voi 
From here, we get that light-like billiard trajectories are 8-periodic in ellipses with 
the ratio of the axes equal to 1 + v2- 
Condition \3.5\) gives: 

y^e |taja|,tany| = {y/2- 1,72 + 1}. 

A few of such trajectories are shown on Figure 



NX •. / 



X .»'' /' 
'• L NX " 



Figure 6. Light- like billiard trajectories with period 8 in the el- 
lipse satisfying a = (3 + 2\/2)b. 
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Example 3.11 (10-perodic light-like trajectories). For n = 10, the condition {5. 



det 



/ B3 Bi B5 Be \ 

Bi B 5 B 6 B 7 

B 5 B e B 7 B s 

\ Bq B 7 B s Bg J 



(a + 6) 20 (5a 2 - lOab + b 2 ){a 2 - lOab + 5b 2 ) 



0. 



(4a6) 22 

From here, we get that light-like billiard trajectories are 10-periodic in ellipses with 
the ratio of the axes equal to either ^1 + — y= or 

From condition \S.5\l , we get: 

a ( ir 2tt 3tt Att 
— G < tan — , tan — , tan — , tan — 
b I 10 10' 10' 10 



Since 



tan- = Wl 



1 



1 



3tt 
tan — 
10 



%/5 V5 + 2V5 tan if 
2 1 1 



V5 ^5-2^5 tan if ' 
both conditions give the same result. 

A few of such trajectories are shown on Figures [7] and [SI 




Figure 7. Light-like billiard trajectories with period 10 in the 
ellipse satisfying a = (1 + 2/\/E)b. 



Example 3.12 (12-perodic light-like trajectories). For n — 12, the condition VS. 3\) 
is: 



det 



( B 3 


Bi 


B 5 


B 6 


B 7 




Bi 


B 5 


B {i 


B- 


B H 




B 5 


Be 


B 7 


B s 


B d 




B 6 


B 7 


Bs 


B<> 


Bio 




\ B 7 


Bs 


B 9 


Bio 


B u 


) 



= 0, 



which can be transformed into: 

{a - 3b)(a - b)('3a -b)(a+ b) 30 (a 2 - Uab + b 2 ) 



(Aab) 22 Vab 



0. 
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Figure 8. Light-like billiard trajectories with period 10 in the 
ellipse satisfying a = (5 + 2y/E)b. 

From here, we get that light-like billiard trajectories are 12-periodic in ellipses with 
the ratio of the axes equal to 2 + 

On the other hand, condition VS. 5}) gives: 

\fl e { tan T^' tan if } ={2- V^' 2 + V3}- 

A few of such trajectories, obtained by a computer simulation, are shown on 
Figure 




Figure 9. Light-like billiard trajectories with period 12 in the 
ellipse satisfying a = (7 + 4V3)6. 

Light-like trajectories in ellipses and rectangular billiards. 

Theorem 3.13. The flow of light-like billiard trajectories within ellipse £ is trajec- 
torially equivalent to the flow of those billiard trajectories within a rectangle whose 
angle with the sides is — . The ratio of the sides of the rectangle is equal to: 

2arc tan 

Remark 3.14. The flow of light-light billiard trajectories within a given oval in 
the Minkowski plane will be trajectorially equivalent to the flow of certain trajecto- 
ries within a rectangle whenever invariant measure m on the oval exists such that 
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m(AB) = m{CD) and m(BC) = m(AD), where A, B, C , D are points on the oval 
where the tangents are light-like. 

Topological properties of elliptical billiards in the Minkowski plane. Here, we are 
going to present the topological description of elliptical billiard in the Minkowski 
plane. In order to do this, we use Fomenko invariants, sec [BF2004] and refer- 
ences therein. The corresponding analysis for the Euclidean case has been done in 
[DR2010] . 

Theorem 3.15. The isoenergy manifold corresponding to the billiard system within 
ellipse £ \3.1}) in the Minkowski plane is represented by the Fomenko graph on 
Figure \1(A 




I 1 1 1 1 

A = +0 \ = a A = 00 X = -b A = -0 

Figure 10. Fomenko graph for elliptical billiard in the Minkowski plane. 

3.2. Relativistic quadrics. The aim of this Section is to present relativistic quad- 
rics, a very recent object, which has been a main tool in our study of billiard dy- 
namics, see jDR2012a) . The point is that geometric quadrics of a confocal pencil 
and their types in pseudo-Euclidean spaces do not satisfy usual properties of con- 
focal quadrics in Euclidean spaces, including those necessary for applications in 
billiard dynamics. For example, we have already mentioned, in the ci-dimcnsional 
Euclidean space, there are d geometric types of quadrics, while in d-dimensional 
pseudo-Euclidean space, there are d + 1 geometric types of quadrics. Thus, we 
first select those important properties of confocal families in the Euclidean spaces 
and axiomatize them as E1-E5. Then we study geometrical types of quadrics in a 
confocal family in the three-dimensional Minkowski space. Next, we analyze tropic 
curves on quadrics in the three-dimensional case and we introduce an important 
notion of discriminant sets T, ± corresponding to a confocal family. The main facts 
about discriminant sets we present in Propositions 13.161 13.181 13.201 13.221 Then, we 
study curved tetrahedra T , which represent singularity sets of E ± and we collect 
related results in Proposition [37TTJ As the next important step, we introduce dec- 
orated Jacobi coordinates for the three-dimensional Minkowski space, and we give 
a detailed description of the colouring into three colours. Each colour corresponds 
to a relativistic type, and we describe decomposition of a geometric quadric of each 
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of the four geometric types into relativistic quadrics. This appears to be a rather 
involved combinatorial-geometric problem, and we solve it by using previous analy- 
sis of discriminant surfaces. We give a complete description of all three relativistic 
types of quadrics. At the end, we generalize definition of decorated Jacobi coordi- 
nates in arbitrary dimensions, and, finally, in Proposition 13.321 we prove properties 
PE1 and PE2, the pseudo-Euclidean analogues of El and E2. 

Confocal quadrics in the three-dimensional Minkowski space and their geometrical 
types. Let us start with the three-dimensional Minkowski space E 2,1 . A general 
confocal family of quadrics in this space is given by: 

(3.6) Q X : + J?- + - J- = 1, A e R, 

a— A o — A c+A 

with a > b > 0, c> 0. 

The family (|3.6j) contains four geometrical types of quadrics: 

• 1-sheeted hyperboloids oriented along z-axis, for A G (-co, — c); 

• ellipsoids, corresponding to A € (— c, b); 

• 1-sheeted hyperboloids oriented along y-axis, for A € (b, a); 

• 2-sheeted hyperboloids, for A 6 (a, +oo) - these hyperboloids are oriented 
along z-axis. 

In addition, there are four degenerated quadrics: Q a , Qb, Q- c , Qoo, that is planes 
x = 0, y = 0, z = 0, and the plane at the infinity respectively. In the coordinate 
planes, we single out the following conies: 

2 2 

y z 

• hyperbola C^ z : ! = 1 in the plane x = 0; 

a — b c + a 

• ellipse Cf z : 1 = 1 in the plane y = 0; 

a — b c + b 

x^ 2 

• ellipse Cf " : V — = 1 in the plane z = 0. 

a + c b + c 

Tropic curves on quadrics in the three-dimensional Minkowski space and discrimi- 
nant sets T,^ . On each quadric, notice the tropic curves - the set of points where 
the induced metrics on the tangent plane is degenerate. 

Since the tangent plane at point (xo, yoy z a) of Q\ is given by the equation: 

xx yy zz Q _ ^ 
a — A 6 — A c + A 
and the induced metric is degenerate if and only if the parallel plane that contains 
the origin is tangent to the light-like cone x 2 + y 2 — z 2 = 0, i.e.: 

2 2 2 
X , Vo Z _ r, 

(a-A) 2 (6-A) 2 (c + A) 2 
we get that the tropic curves on Q\ are the intersection of the quadric with the 



cone: 

,2 



y =o, 



(a - A) 2 (b - A) 2 (c + A) 2 
see |KT2Q09j . 

Now, consider the set of the tropic curves on all quadrics of the family (|3.6p . 
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Proposition 3.16. The union of the tropic curves on all quadrics of LS.6\) is a 

union of two ruled surfaces S + and £~ which can be parametrically represented as: 



a — X b — A I cos 2 t sin 2 1 

■.cost, y = — ■== sin t, z — {c + X)] 



\Ja + c ' \Jb + c \a + c 6 + c' 



a — X b — X I cos 2 t sin 2 t 

cost, y = —r^^= smt, z = — (c + Ah 



^/a + c Vb~+c ya + c b + c 

with X e R, t e [0,2vr). 

The intersection of these two surfaces is an ellipse in the xy-plane: 



,2 



£+n£- : + -2_ = 1, z = 0. 

a + c o + c 

TTie too surfaces S + , £~ are developable as embedded into Euclidean space. More- 
over, their generatrices are all light-like. 

Surfaces S + and £~ from Proposition 13. 161 are represented on Figure [TO 



Figure 1 1 . The union of all tropic curves of a confocal family. 



In |Peil999| , a definition of generalization of Gauss map to surfaces in the three- 
dimensional Minkowski space is suggested. Namely, the pseudo vector product is 
introduced as: 



(x 2 y3-x 3 y2, x 3 yi-xiy 3 , -(x 1 y 2 -x 2 yi)) = £ 2 ,i(xxy). 





Xl 


X2 


%3 


xAy = 


yi 


V2 


2/3 




ei 


e 2 


-e3 



It is easy to check that (x A y, x) 2 1 = (xAyj) 21 = 0. 
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Then, for surface S : U — > E 2 ' 1 , with U C R 2 , the Minkowski Gauss map is 
defined as: 



where P : R 3 \ {(0, 0, 0)} — > RP 2 is the usual projectivization. 

Since r\ A r t is light-like for all A and t, the Minkowski Gauss map of surfaces 
S ± is singular at all points. 

The pseudo-normal vectors to are all light-like, thus these surfaces are light- 
like developable, as defined in [CI2010] , There, a classification of such surfaces is 
given - each is one part-by-part contained in the following: 

• a light-like plane; 

• a light-like cone; 

• a tangent surface of a light-like curve. 

Since E + and E - are contained neither in a plane nor in a cone, we expect that 
they will be tangent surfaces of some light-like curve, which is going to be shown 
in the sequel, see Corollary 13.211 later in this section. 

On each of the surfaces S + , we can notice that tropic lines corresponding 
to 1-sheeted hyperboloids oriented along y-axies form one curved tetrahedron, see 
Figure 1111 Denote the tetrahedra by T + and T~ respectively: they are symmetric 
with respect to the xy-plane. On Figure [12] tetrahedron T + C S + is shown. 



Figure 12. Curved thetrahedron T + : the union of all tropic cur- 
ves on E + corresponding to A € (b, a). 



Let us summarize the properties of these tetrahedra. 

Proposition 3.17. Consider the subset T + of E + determined by the condition 
A 6 [b, a] . This set is a curved tetrahedron, with the following properties: 
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• its verteces are: 



r,|4=L o, 4+LV vJ JULu. ; ' ' 



y/a + c ' yja + cj ' \ yja + c ' yja + c 
a — b a + c \ / a — b a + c 

y 3 o, -==, -== , v 4 o, 



y/b + c VbTcJ ' V ' y/b + c Vb + 

• f/ie shorter arcs of conies C^ z and C^ z determined by V\, V2 and V3, V4 
respectively are two edges of the tetrahedron; 

• those two edges represent self-intersection o/S + ; 

• other four edges are determined by the relation: 

(3.7) -a-6 + 2A+(a-&)cos2< = 0, 

• those four edges are cuspidal edges o/S + ; 

• thus, at each vertex of the tetrahedron, a swallowtail singularity of S + oc- 
curs. 

It can be proved that the tropic curves of the quadric Q\ represent exactly the 
locus of points [x, y, z) where equation 

2 ? ? 

x y z 

y ' a-X b-X c+X 

has Ao as a multiple root. 

Proposition 3.18. A tangent line to the tropic curve of a non-degenarate quadric 
of the family h3. 6\) is always space-like, except on a 1-sheeted hyperboloid oriented 
along y-axis. 

Tangent lines of a tropic on l-sheeted hyperboloids oriented along y-axis are light- 
like exactely at four points, while at other points of the tropic curve, the tangents 
are space-like. 

Moreover, a tangent line to the tropic of a quadric from \3. 6\) belongs to the 
quadric if and only if it is light-like. 



Remark 3.19. In other words, the only quadrics of the family i3. 6}) that may 
contain a tangent to its tropic curve are l-sheeted hyperboloids oriented along y- 
axis, and those tangents are always light-like. The tropic curves and their light-like 
tangents on such an hyperboloid are shown on Figure \73l 

Notice that equations obtained in the proof of Proposion 13.181 are equivalent to 
equation (|3.7p of Proposition ^. 171 which leads to the following: 

Proposition 3.20. Each generatrix of E + and X~ is contained in one l-sheeted 
hyperboloid oriented along y-axis from i3.6\) . Moreover, such a generatrix is touch- 
ing at the same point one of the tropic curves of the hyperboloid and one of the 
cusp-like edges of the corresponding curved tetrahedron. 

Corollary 3.21. Surfaces E + and £~ are tangent surfaces of the cuspidal edges 
of thetrahedra T + and T~ respectively. 

In next propositions, we give further analysis the light-like tangents to the tropic 
curves on an l-sheeted hyperboloid oriented along y-axis. 
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Figure 13. The tropic curves and its light-like tangents on a hyperboloid. 

Proposition 3.22. For a fixed Ao € (b,a), consider a hyperboloid Q\ from H3. b}) 
and an arbitrary point (x,y,z) on Q\ . Equation \3.8)) has, along with Ao, two 
other roots in C : denote them by Xi and A2 . Then Ai =A2 if and only if (x,y,z) is 
placed on a light-like tangent to a tropic curve of Q\ . 

Proof. Follows from the fact that the light-like tangents are contained in the E + U 
£~, see Propositions 13.161 and l3.2l)l □ 

Proposition 3.23. Two light-like lines on a one-sheeted hyperboloid oriented along 
y-axis from i3. 6}) are either skew or intersect each other on a degenerate quadric 
from fttyl 

Proof. Follows from the fact that the hyperboloid is symmetric with respect to the 
coordinate planes. □ 

Lemma 3.24. Consider a non-degenerate quadric Q\ , which is not a hyperboloid 
oriented along y-axis, i.e. Ao $jL [b,a] U {— c}. Then each point of Q\ which is not 
on one of the tropic curves is contained in two additional distinct quadrics from the 
family 13.6}) . 

Consider two points A, B of Q\ , which are placed in the same connected compo- 
nent bounded by the tropic curves, and denote by X' A , A^ and \' B , Xg the solutions, 
different than Ao, of equation \3. ffj) corresponding to A and B respectively. Then, 
if Xq is smaller (resp. bigger, between) than X' A , X A , it is also smaller (resp. bigger, 
between) than X' B , Xg. 

Lemma 3.25. Let Q\ be a hyperboloid oriented along y-axis, Xq € (b,a), and A, 
B two points of Q\ a , which are placed in the same connected component bounded 
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by the tropic curves and light-like tangents. Then, if A is contained in two more 
quadrics from the family 13. 6\) . the same is true for B. 

In this case, denote by X' A , X A and X' B , X B the real solutions, different than A , 
of equation A3.8\) corresponding to A and B respectively. Then, if X is smaller 
(resp. bigger, between) than X' A , X A , it is also smaller (resp. bigger, between) than 
X' B , X B . 

On the other hand, if A is not contained in any other quadric from \3.b}) . then 
the same is true for all points of its connected component. 

Proof. The proof of both Lemmae 13.241 and 13.251 follows from the fact that the 
solutions of (|3.8|) are continuously changed through the space and that two of the 
solutions coincide exactly on tropic curves and their light-like tangents. □ 

Generalized Jacobi coordinates and relativistic quadrics in the three-dimensional 
Minkowski space. 

Definition 3.26. Generalized Jacobi coordinates of point (x,y,z) in the three-di- 
mensional Minkowski space E 2 ' 1 is the unordered triplet of solutions of equation 

Note that any of the following cases may take place: 

• generalized Jacobi coordinates are real and different; 

• only one generalized Jacobi coordinate is real; 

• generalized Jacobi coordinates are real, but two of them coincide; 

• all three generalized Jacobi coordinates are equal. 

Lemmae 13.241 and 13.251 will help us to define relativistic types of quadrics in 
the 3-dimensional Minkowski space. Consider connected components of quadrics 
from (|3.6|) bounded by tropic curves and, for 1-sheeted hyperboloids oriented along 
y-axis, their light-light tangent lines. Each connected component will represent a 
relativistic quadric. 

Definition 3.27. A component of quadric Q\ a is of relativistic type E if, at each 
of its points, Xq is smaller than the other two generalized Jacobi coordinates. 

A component of quadric Q\ n is of relativistic type H 1 if, at each of its points, 
Ao is between the other two generalized Jacobi coordinates. 

A component of quadric Q\ is of relativistic type H 2 if, at each of its points, 
Ao is bigger than the other two generalized Jacobi coordinates. 

A component of quadric Q\ a is of relativistic type if, at each of its points, Ao 
is the only real generalized Jacobi coordinate. 

Lemmae 13.241 and 13.251 guarantee that types of relativistic quadrics are well- 
defined, i.e. that to each such a quadric a unique type E, H 1 , H 2 , or can be 
assigned. 

Definition 3.28. Suppose (x,y,z) is a point of the three-dimensional Minkowski 
space E 2,1 where equation \ 3.8\ has real and different solutions. Decorated Jacobi 
coordinates of that point is the ordered triplet of pairs: 

(£,Ai), (H\X 2 ), (H 2 ,X 3 ), 

of generalized Jacobi coordinates and the corresponding types of relativistic quadrics. 
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Now, we are going to analyze the arrangement of the relativistic quadrics. Let 
us start with their intersections with the coordinate planes. 

Intersection with the xy-plane. In the xy-plane, the Minkowski metrics is reduced 
to the Euclidean one. The family (|3.6|) is intersecting this plane by the following 
family of confocal conies: 



(3.9) 

see Figure ITU 



y 



A b-X 



= 1, 




Figure 14. Intersection of relativistic quadrics with the xy-plane. 



We conclude that the xy-plane is divided by ellipse into two relativistic 
quadrics: 

• the region within is a relativistic quadric of -E-type; 

• the region outside this ellipse is of if 1 -type. 

Moreover, the types of relativic quadrics intersecting the xy-plane are: 

• the components of ellipsoids are of H 1 -type; 

• the components of 1-sheeted hyperboloids oriented along y-axis of H 2 -type; 

• the components of 1-sheeted hyperboloids oriented along z-axis of -E-type. 

On Figure type E quadrics are coloured in dark gray, type H 1 medium gray, 
and type H 2 light gray. The same colouring rule is applied on Figures [T51IT71 
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Intersection with the xz-plane. In the xz-plane, the reduced metrics is the Min- 
kowski one. The intersection of family (|3.6j) with this plane is the following family 
of confocal conies: 



(3.10) 

see Figure IT5l 



a — A c + A 



1. 




Figure 15. Intersection of relativistic quadrics with the xz-plane. 



The plane is divided by ellipse C% z and the four joint tangents of (|3.10[) into 13 
parts: 

• the part within C£ z is a relativistic quadric of i/ 1 -type; 

• four parts placed outside of C£ z that have non-empty intersection with the 
x-axis are of i? 2 -type; 

• four parts placed outside of C£ z that have non-empty intersection with the 
z-axis are of _B-type; 

• the four remaining parts are of 0-type and no quadric from the family (|3 - 6[) , 
except the degenerated Qb, is passing through any of their points. 

Intersection with the yz-plane. As in the previous case, in the yz-plane, the reduced 
metrics is the Minkowski one. The intersection of family Q3.6P with this plane is 
the following family of confocal conies: 



(3.11) 

see Figure [TBI 



r>v z 



b - A c + A 



= 1, 
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Figure 16. Intersection of relativistic quadrics with the j/z-plane. 

The plane is divided by hyperbola C^ z and joint tangents of (|3 . 11 1) into 15 parts: 

• the two convex parts determined by C% z are relativistic quadric of H 1 -type; 

• five parts placed outside of C^ z that have non-empty intersection with the 
coordinate axes are of ff 2 -type; 

• four parts, each one placed between C^ z and one of the joint tangents of 
flnHU) are of S-type; 

• through points of the four remaining parts no quadric from the family (|3.6[) , 
except the degenerated Q a , is passing. 

Intersection of relativistic quadrics with the coordinate planes is shown in Figure 

Let us notice that from the above analysis, using Lemmae 13.241 and 13.251 we can 
determine the type of each relativistic quadric with a non-empty intersection with 
some of the coordinate hyper-planes. 

1-sheeted hyperboloids oriented along z-axis: A € (— oo, — c). Such a hyperboloid is 
divided by its tropic curves into three connected components - two of them are 
unbounded and mutually symmetric with respect to the xy-plane, while the third 
one is the bounded annulus placed between them. The two symmetric ones are of 
i? 1 -type, while the third one is of -E-type. 

Ellipsoids: A £ (— c, 6). An ellipsoid is divided by the tropic curves into three 
bounded connected components - two of them are mutually symmetric with respect 
to the ccy-plane, while the third one is the annulus placed between them. In this 
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FIGURE 17. Intersection of relativistic quadrics with coordinate planes. 

case, the symmetric components represent relativistic quadrics of -E-type. The 
annulus is of 7J 1 -type. 

1-sheeted hyperboloids oriented along y-axis: A € (b, a). The decomposition of those 
hyperboloids into relativistic quadrics is more complicated and interesting than for 
the other types of quadrics from (|3.6j) . By its two tropic curves and their eight light- 
like tangent lines, such a hyperboloid is divided into 28 connected components: 

• two bounded components placed inside the tropic curves are of H -type; 

• four bounded components placed between the tropic curves and light-like 
tangents, such that they have non-empty intersections with xz-plane are of 
-ff 2 -type; 
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• four bounded components placed between the tropic curves and light-like 
tangents, such that they have non-empty intersections with yz-plane are of 
S-type; 

• two bounded components, each limited by four light-like tangents, are of 
# 2 -type; 

• four unbounded components, each limited by two light-like tangents, such 
that they have non-empty intersections with the xy-plane, are of if 2 -type; 

• four unbounded components, each limited by two light-like tangents, such 
that they have non-empty intersections with the j/z-plane, are of -E-type; 

• eight unbounded components, each limited by four light-like tangents, are 
sets of points not contained in any other quadric from (|3.6[) . 

2-sheeted hyperboloids: X € (a, +00). Such a hyperboloid is by its tropic curves 
divided into four connected components: two bounded ones are of i/ 2 -type, while 
the two unbounded are of if -type. 

Decorated Jacobi coordinates and relativistic quadrics in d- dimensional pseudo- 
Euclidean space. Now we are going to introduce relativistic quadrics and their types 
in confocal family (|2.4|) in the d-dimensional pseudo-Euclidean space E fe,/ . 

Definition 3.29. Generalized Jacobi coordinates of point x in the d- dimensional 
pseudo-Euclidean space E fc '' is the unordered d-tuple of solutions X of equation: 

(3.12) -J^. + ... + _g^+ +... + _g;L. = i. 

01 - A a k - X a k+1 + X aj, + X 

As already mentioned in Section 12.31 this equation has either d or d — 2 real 
solutions. Besides, some of the solutions may be multiple. 

The set E^ of points x in R d where equation (I3.12[) has multiple solutions is an 
algebraic hyper-suface. E^ divides each quadric from (|2.4[) into several connected 
components. We call these components relativistic quadrics. 

Since the generalized Jacobi coordinates depend continuosly on x, the following 
definition can be made: 

Definition 3.30. We say that a relativistic quadric placed on Q\ g is of type E if, at 
each of its points, Xq is smaller than the other d—1 generalized Jacobi coordinates. 

We say that a relativistic quadric placed on Q\ is of type H l (1 < % < d — 1) if, 
at each of its points, Xq is greater than other i generalized Jacobi coordinates, and 
smaller than d — i — 1 of them. 

We say that a relativistic quadric placed on Q\ is of type l (0 < i < d — 2) if, 
at each of its points, Xq is greater than other i real generalized Jacobi coordinates, 
and smaller than d — i — 2 of them. 

It would be interesting to analyze properties of the discriminant manifold E^, as 
well as the combinatorial structure of the arrangement of relativistic quadrics, as 
it is done for d = 3. Remark that this description would have [d/2] substantially 
different cases in each dimension, depending on choice of k and I. 

Definition 3.31. Suppose (xi, . . . , X4) is a point of the d-dimensional Minkowski 
space E fc '' where equation \3.12\) has real and different solutions. Decorated Jacobi 
coordinates of that point is the ordered d-tuplet of pairs: 

(E,Xx), (H\X 2 ), (H d -\X d ), 
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of generalized Jacobi coordinates and the corresponding types of relativistic quadrics. 

Since we will consider billiard system within ellipsoids in the pseudo-Euclidean 
space, it is of interest to analyze behaviour of decorated Jacobi coordinates inside 
an ellipsoid. 

Proposition 3.32. Let £ be ellipsoid in E fe,/ given by S2.3\) . We have: 

PE1 each point inside £ is the intersection of exactly d quadrics from ^2.4\ l; 

moreover, all these quadrics are of different relativistic types; 

PE2 the types of these quadrics are E, H , . . . , - each type corresponds 

to one of the disjoint intervals of the parameter A: 

(-a d ,-a d ^i), (—a,d-i, —0,4-2), • (-Ofc+i,0), (0,Ofe), (ofe,Ofc_i), (02,01). 

Proof. The function given by the left-hand side of (|3.12[) is continous and strictly 
monotonous in each interval (-a<j, —aa-i), (-0^-1,-0^-2), (-o*+2, -0^+1), 
(ofc,Ofe_i), (02,01) with infinite values at their cndpoints. Thus, equation 
(|3.12[) has one solution in each of them. On the other hand, in (— a^, +1 ,afe), the 
function is tending to +00 at the endpoints, and has only one extreme value - 
the minimum. Since the value of the function for A = is less than 1 for a point 
inside £, it follows that equation (13.121) will have two solutions in (— a/c+i, a^) - one 
positive and one negative. □ 

3.3. Billiards within quadrics and their periodic trajectories. In this sec- 
tion, we are going to derive first further properties of ellipsoidal billiards in the 
pseudo-Euclidean spaces. First we find in Theorem 13.331 a simple and effective cri- 
terion for determining the type of a billiard trajectory, knowing its caustics. Then 
we derive properties PE3-PE5 in Propostion 13.341 After that, we prove the gen- 
eralization of Poncelet theorem for ellipsoidal billiards in pseudo-Euclidean spaces 
and derive the corresponding Cayley-type conditions, giving a complete analytical 
description of periodic billiard trajectories in arbitrary dimension. These results 
are contained in Theorems 13.351 and 13.361 

Ellipsoidal billiards. 

Ellipsoidal billiard. Billiard motion within an ellipsoid in the pseudo-Euclidean 
space is a motion which is uniformly straightforward inside the ellipsoid, and obeys 
the reflection law on the boundary. Further, we will consider billiard motion within 
ellipsoid £, given by equation (|2.3|) in ~E k ' L . The family of quadrics confocal with £ 
is (12T)) . 

Since functions Fi given by (|2.8|) are integrals of the billiard motion (see Mosl980, 
IAudl994|iKT2009| ) . we have that for each zero A of the equation (1277]) . the corre- 
sponding quadric Q\ is a caustic of the billiard motion, i.e. it is tangent to each 
segment of the billiard trajectory passing through the point x with the velocity 
vector v. 

Note that, according to Theorem 12.181 for a point placed inside £, there are 
d real solutions of equation (13. 12[) . In other words, there are d quadrics from the 
family (|2.4p containing such a point, although some of them may be multiple. Also, 
by Proposition 12.171 and Theorem [233 a billiard trajectory within an ellipsoid will 
always have d — 1 caustics. 
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According to Remark I2.16[ all segments of a billiard trajectory within £ will be 
of the same type. Now, we can apply the reasoning from Section 12.31 to billiard 
trajectories: 

Theorem 3.33. In the d- dimensional pseudo- Euclidean space E fc, ' ; consider a bil- 
liard trajectory within ellipsoid £ = Q®, and let quadrics Q ai , ■ Qa d _i from 
the family \2.J$ be its caustics. Then all billiard trajectories within £ sharing the 
same caustics are of the same type: space-like, time-like, or light-like, as the initial 
trajectory. Moreover, the type is determined as follows: 

• if oo G {a\, . . . ,ctd-i}, the trajectories are light-like; 

• if (— 1) • a.\ ■ . . . ■ ctd-i > 0, the trajectories are space-like; 

• if (—1)' • Q.\ ■ ■ ■ ■ ■ a<j_i < 0, the trajectories are time-like. 

Proof. Since values of functions Fi given by (|2.8[) are preserved by the billiard 
reflection and 

d 

^2Fi(x,v) = (v,v) k j, 

i=l 

the type of the billiard trajectory depends on the sign of the sum 2i=i F%{ x i v )- 
From the equivalence of relations (|2.7[) and (12.9[) . it follows that the sum depends 
only of the roots of V, i.e. of parameters a±, . . . , a^-i of the caustics. 

Notice that the product a.\ ■ . . . ■ otd-i is changed continuously on the variety of 
lines in E fc: ' that intersect £, with infinite singularities at light-like lines. Besides, 
the subvariety of light-like lines divides the variety of all lines into subsets of space- 
like and time-like ones. When passing through light-like lines, one of parameters cti 
will pass through the infinity from positive to the negative part of the reals or vice 
versa; thus, a change of sign of the product occurs simultaneously with a change of 
the type of line. 

Now, take otj = —ak+j for 1 < j < I, and notice that all lines placed in the 
/c-dimensional coordinate subspace E fe x Z will have the corresponding degenerate 
caustics. The reduced metrics is Euclidean in this subspace, thus such lines are 
space-like. Since cti, . . . , otk are positive for those lines of E* x Z that intersect £, 
the statement is proved. □ 

Let us note that, in general, for the fixed d— 1 quadrics from the confocal family, 
there can be found joint tangents of different types, which makes Theorem [3T33] in a 
way unexpected. However, it turns out that, with fixed caustics, only lines having 
one type may have intersection with a given ellipsoid — and only these lines give 
rise to billiard trajectories. 

Next, we are going to investigate the behaviour of decorated Jacobi coordinates 
along ellipsoidal billiard trajectories. 

Proposition 3.34. LetT be a trajectory of the billiard within ellipsoid £ in pseudo- 
Euclidean space E fe,z . Denote by ai, . . . , a^-i the parameters of the caustics from 
the confocal family \2.J$ ofT, and take b\, . . . , b p , c\, . . . , c q as in Theorem \2.18[ 
Then we have: 

PE3 along T , each generalized Jacobi coordinate takes values in exactly one of 
the segments: 

[c 2 /-i,c 2 ;- 2 ], [c 2 ,ci], [ci,0], [0,6i], [b 2 ,b 3 }, [b 2 k-2, hk-i}; 



VLADIMIR DRAGOVIC AND MILENA RADNOVIC 



PE4 along T , each generalized Jacobi coordinate can achieve local minima and 
maxima only at touching points with corresponding caustics, intersection 
points with corresponding coordinate hyper-planes, and at reflection points; 

PE5 values of generalized Jacobi coordinates at critical points are 0, b\, 

&2fc-i; ci, . . . , C2/-1; between the critical points, the coordinates are changed 
monotonously. 



Proof. Property PE3 follows from Theorem 12.181 Along each line, the generalized 
Jacobi coordinates are changed continuously. Moreover, they are monotonous at 
all points where the line has a transversal intersection with a non-degenerate quad- 
ric. Thus, critical points on a line are exactly touching points with corresponding 
caustics and intersection points with corresponding coordinate hyper-planes. 

Note that reflection points of T are also points of transversal intersection with 
all quadrics containing those points, except with £ . Thus, at such points, will be 
a critical value of the corresponding generalized Jacobi coordinate, and all other 
coordinates are monotonous. This proves PE4 and PES. □ 



The properties we obtained are pseudo-Euclidean analogs of properties E3-E5, 
which are true for ellipsoidal billiards in Euclidean spaces. 



Analytic conditions for periodic trajectories. Now, we are going to derive the cor- 
responding analytic conditions of Cayley's type for periodic trajectories of the el- 
lipsoidal billiard in the pseudo-Euclidean space, and therefore to obtain the gener- 
alization of the Poncelet theorem to pseudo-Euclidean spaces. 

Theorem 3.35 (Generalized Cayley-type conditions). In the pseudo-Euclidean 
space E k ' 1 (k + I = d), consider a billiard trajectory T within ellipsoid £ given 
by equation \2.3\) . Let Q ai , ■ ■ ■ , Qa d _i from confocal family {2.1$ be caustics ofT- 
Then T is periodic with period n if and only if the following condition is satisfied: 

• for n = 2m: 
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for n = 2m + 1 : 
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Here, (-Bj), (Ci), (-Dj) are coefficients in the Taylor expansions around X = of the 

/(A) 



functions /(A) = \J (a± — A) 

— - respectively. 
C\ — A 

Proof. Denote: 

^i(A) = (ax 



(a d -i - A) • (at - e 1 X) 



(a d - e d X), 



h - A' 



A) • . . . • (ad-i " A) • (ai - eiA) 



(a d - e d X). 



Following Jacobi [Jacl884 , along a given billiard trajectory, we consider the inte- 
grals: 



(3.13) 



d 

E 

s=l 



dX, 



E 



AadAs 



d 

E 

s=l 



x d s - 2 dx s 



By PE3 of Proposition 13.341 we may suppose that: 

Ai G [0,6i], Aj G [62t-2,i>2i-i] for 2 < i < fc; 
Afc+i G [ci,0], Xk+j G [c2j-i,C2i-2] for 2 < j <l. 

Along a billiard trajectory, by PE4 and PE5 of Proposition 13.341 each X s will pass 
through the corresponding interval monotonously from one endpoint to another 
and vice versa alternately. Notice also that values b\, 62/t-i, c%, ou—x 
correspond to the Weierstrass points of hyper-elliptic curve: 

(3.14) M 2 =7MA). 

Thus, calculating integrals (|3.13p . we get that the billiard trajectory is closed after 
n reflections if and only if, for some Tlx, n% such that ni + n-i = n: 

nA(Po) = n 1 A(P bl )+n 2 A(P Cl ), 



on the Jacobian of curve (|3.14|) . Here, A is the Abel- Jacobi map, and Pt is a point 
on the curve corresponding to A = t. Further, in the same manner as in [DR1998b , 
we obtain the conditions as stated in the theorem. □ 



As an immediate consequence, we get: 

Theorem 3.36 (Generalized Poncelet theorem). In pseudo-Euclidean space E fe -' 
(k + I = d), consider a billiard trajectory T within ellipsoid £ . 
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If T is periodic and become closed after n reflections on the ellipsoid, then any 
other trajectory within £ having the same caustics as T is also periodic with period 
n. 

Remark 3.37. The generalization of the Full Poncelet theorem from [CCS1993] to 

pseudo-Euclidean spaces is obtained in |WFS + 2009 . However, only space-like and 



time-like trajectories were discussed there. 

A Poncelet-type theorem for light-like geodesies on the ellipsoid in the three- 
dimensional Minkowski space is proved in [GKT2007] . 

Remark 3.38. Theorems \3.35\ and \3.36\ will also hold in symmetric and degenerated 
cases, that is when some of the parameters SiOi, ctj concide, or in the case of light- 
like trajectories, when oo € {ay \ 1 < j < d — 1}. In such cases, we need to apply 
the desingularisation of the corresponding curve, as explained in detail in our works 
[DR2006bllDR2008| . 

When we consider light-like trajectories, then the factor containing the infinite 
parameter is ommited from polynomial V\ . 



Example 3.39. Let us find all ^-periodic trajectories within ellipse £ given by \3. 1\) 
in the Minkowski plane, i.e. all conies C a from the confocal family \3.2\l correspond- 
ing to such trajectories. 

By Theorem \3.35\ the condition is B 3 — 0, with 



y/(a- X)(b + X){a- A) = B + Si A + B 2 X 2 + B 3 X 3 + . 
being the Taylor expansion around X = 0. Since 

(—ab — aa + ba)(—ab + act + ba)(ab + aa + ba) 
3 = 16(ak*) 5/2 ' 

we obtain the following solutions: 

ab ab ab 

ai = -. , a 2 = — — r, «3 = 



b — a' a + 6' ' a + b 

Since a\ $ (—b,a) and 0:2,0:3 € (—6, a), conic C ai is a hyperbola, while C a2 , C a3 
are ellipses. 

4. INTEGRABLE LINE CONGRUENCES AND DOUBLE REFLECTION NETS 

As a modern scientific discipline, the Discrete Differential Geometry emerged 
quite recently (see BS2008 ), within a study of the lattice geometry. So-called 
intcgrability conditions for quad-graphs have a fundamental role there. On the 
other hand, the geodesies on an ellipsoid are one of the most important and exciting 
examples of the classical Differential Geometry. Billiard systems within quadrics are 
known to be seen as natural discretizations of the systems of geodesies on ellipsoids. 
In the sequel we are going to present them as a part of the building blocks of the 
foundations of the Discrete Differential Geometry. 

The main elements of the systems on quad-graphs are the equations of the form 
Q(x, x%, a?2) £12) = on quadrilaterals, where Q is a multiafhne polynomial, that 
is a polynomial of degree one in each argument. Such equations are called quad- 
equations. The field variables Xi are assigned to four vertices of a quadrilateral as in 
Figure [I8l The quad equation can be solved for each variable, and the solution is a 
rational function of the other three variables. Following ABS2009). we consider the 
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Figure 18. Quad-equation Q(xi,X2,X3,X4,) = 0. 
idea of integrability as consistency, see Figure [T9j We assign six quad-equations to 
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X\ 



Figure 19. 3D-consistency. 

the faces of a coordinate cube. The system is said to be ZD -consistent if the three 
values for X123 obtained from the equations on right, back, and top faces coincide 
for arbitrary initial data. 

We will be interested here in a geometric version of the integrable quad graphs, 
with lines in P d playing a role of the vertex fields. We will denote by C d the 
Grassmannian Gr(2,d + 1) of two-dimensional vector subspaces of the (d + 1)- 
dimensional vector space, d > 2. 

4.1. Billiard algebra and quad-graphs. This section is devoted to a quad-graph 
interpretation from |DR2012"b] of some results obtained previously using billiard 
algebra. 

Let us start from a theorem on confocal families of quadrics from DR2008J : 

Theorem 4.1 (Six-pointed star theorem). Let J 7 be a family of confocal quadrics in 
P 3 . There exist configurations consisting of twelve planes in P 3 with the following 
properties: 

• The planes may be organized in eight triplets, such that each plane in a 
triplet is tangent to a different quadric from T and the three touching points 
are collinear. Every plane in the configuration is a member of two triplets. 

• The planes may be organized in six quadruplets, such that the planes in 
each quadruplet belong to a pencil and are tangent to two different quadrics 
from J-. Every plane in the configuration is a member of two quadruplets. 

Moreover, such a configuration is determined by three planes tangent to three dif- 
ferent quadrics from T , with collinear touching points. 



40 



VLADIMIR DRAGOVIC AND MILENA RADNOVIC 



Such a configuration of planes in the dual space P 3 * is shown in Figure [20] each 
plane corresponds to a vertex of the polygonal line. 



X1X12 




X3X23 



Figure 20. A configuration of planes from Theorem th:zvezda 

To understand the notation used in Figure [20] let us recall the construction 
leading to configurations from Theorem 14. II Take Qi, Q2, Q3 to be quadrics from 
J 7 , and a, /3, 7 respectively their tangent planes such that the touching points A, B, 
C are collinear. Denote by x the line containing these three points, and by x\, X2, 
X3 the lines obtained from x by reflections off Qi, Q2, Q3 at A, B, C respectively. 

Now, as in Proposition 12.141 determine lines X12, X13, X23, X123 such that they 
respectively complete triplets {x, x\, X2}, {x,xi,x^}, {x, X2, £3}, {x^, £13, X23} to 
double reflection configurations H 

Notice the following objects in Figure [2(J] 

twelve verteces: to each vertex, a plane tangent to one of the three quad- 
rics Qi, Q2, Q3 and a pair of lines are assigned — the lines of any pair 
are reflected to each other off the quadric at the touching point with the 
assigned plane; 

eight triangles: in any triangle, the planes assigned to the verteces are 
touching the corresponding quadrics at three collinear points — thus to 
each triangle, the line containing these points is naturaly assigned; 

six edges: each edge containes four verteces — four planes assigned to these 
verteces are in the same pencil; thus a double reflection configuration cor- 
responds to each edge. 

Now, we are ready to prove the 3-D-consistency of the quad-relation introduced 
via double reflection configurations. 



3 Let us note that in DR2008 , the lines x, x\, X2, £3, £12, £'13, £'23, 2:123 were respectively 
denoted by 0, p, q, s, —x, pi, qi, x + s, where the addition is defined in the billiard algebra 
introduced in that paper. 
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Theorem 4.2. Let x, x\, x%, X3 be lines in the projective space, such that X\, X2, 
X3 are obtained from x by reflections off confocal quadrics Q\, Q2, Q3 respectively. 
Introduce lines xyi, x\z , X23, £123 such that the following quadruplets are double 
reflection configurations: 

{x,Xi 7 Xi2,X 2 }, {x,Xi,Xi3,Xz}, {x,X2,X23,Xz}, {xi , X 12 , £123 , Zl3 } ■ 

Then the following quadruplets are also double reflection configurations: 

{x 2 , X Y 2,X\23, 2:23}, {^3, ^13, &123, ^23}- 

Proof. Let us remark that the configuration described in Theorem 14.11 has obviously 
a combinatorial structure of a cube, with planes corresponding to the edges of the 
cube. In this way, lines x, x±, X2, £3, £12, £13, £23, X123 will correspond to the 
vertices of the cube as shown in Figure fT9l A pair of lines is represented by endpoints 
of an edge if they reflect to each other off the plane joined to this edge. Faces of the 
cube represent double reflection configurations. Notice also that planes joined to 
parallel edges of the cube are tangent to the same quadric. The statement follows 
from Theorem 14.11 and the construction given after, see Figure [20l □ 

4.2. Double reflection nets. Assume a family of confocal quadrics is given in 
P d . Notice that, by the Chasles theorem [Chal827 , every line in P d touches d — 1 
quadrics from the family. 

Moreover, by Corollary 12. 10[ these d — 1 quadrics are preserved by the billiard 
reflection. Confocal quadrics touched by a line are called caustics of this line, or 
consequently, caustics of the billiard trajectory that contains the line. 

Now, fix d — 1 quadrics from the pencil an take A C C d to be the set of all lines 
touching these d — 1 quadrics. 

Definition 4.3. A double reflection net is a map 

(4.1) p : T n ^A, 

such that there exist m quadrics Qi, . . . , Q m from the confocal pencil, satisfying 
the following conditions: 

(1) sequence {^p{iiq + iej)}; e z represents a billiard trajectory within Qj , for 
each j € {1, ... , m} and n € Z m ; 

(2) lines ip(no), (p(no + ei), (p(no + ej) 7 ip(iio + ej + e^) form a double reflection 
configuration, for all i, j G {1, . . . , to}, i 7^ j and no € Z m . 

In other words, for each edge in Z m of direction e,, the lines corresponding to 
its vertices intersect at Qi, while the four tangent planes at the intersection points, 
associated to an elementary quadrilateral, belong to a pencil. 

In the following subsections, we describe some examples of double reflection nets. 
After that, we construct F-transformations of double reflection nets and conclude 
this section by establishing connection with the Grassmanian Darboux nets from 
[AHS2009| . 

Example of a double reflection net in the Minkowski space. Consider three-dimensional 
Minkowski space E 2 ' 1 . In this space, let a general confocal family is given by (|3.6|) . 

Fix Ao € (6, a), and consider hyperboloid Q\ . Denote by 01, a 2 , 03, a 4 , b\, b 2 , 
63, 64 the light-like generatrices of Q\ , in the following way (see Figure [211 : 
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FIGURE 21. The tropic curves of Q\ and their light-like tangents. 



• lines en belong to one, and hi to the other family of generatrices of Q\ ; 
that is, di and a,j are always skew for i ^ j, while dj and bj are coplanar 
for all 

• oi, a%, 63, 64 are tangent to the tropic curve contained in the half-space 
z > 0, while a^, a^, 61, 62 are touching the other tropic curve; 

• dj is parallel to hi for each i; 

• pairs (01,62), (0.2, bi), (03,64), (04,63) have intersection points in the xy- 
plane; 

• pairs (04,63), (02,64), (03,61), (04,62) have intersection points in the xz- 
plane; 

• pairs (01,64), (02,63), (03,62), (04,61) have intersection points in the yz- 
plane. 

Take A to be the set of all generatrices of hyperboloid Q\ , i.e. the set of all 
lines having Q\ as a double caustic. In particular, A contains all lines Oj, 6j. 
It is possible to define a map 



<Pm ■ Z^^A, 
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such that the image of ifiM is the set {ai, 02, 03, 04, b\, 62, 03, 64} and for each neZ 4 
lines ^m(ii + ei), </2&r(n + e 2 ), fM (n + e 3 ), <^jw(n + 64) are obtained from <£m(ii) 
by reflection off Q OJ Q b , Q_ c , Qoo respectively. 

More precisely, (pM will be periodic with period 2 in each coordinate and: 

(0,0, 0,0) =<pu(l, 1,1,1) =oi, ^ M (l,l,0,0) = ^ M (0,0,l,l) = a 2 , 
¥»*r(l, 0, 1, 0) = ^m(0, 1,0, 1) = a 3 , ^ M (0, 1, 1, 0) = <p M (l, 0, 0, f ) = a 4 , 
^m(0, 0, 0, f ) = </? M (l, 1, 1,0) = h, <p M {l, 1, 0, 1) = ^m(0, 0, 1,0) = b 2 , 
^(1,0,1,1) = <p M (0, 1,0,0) = a 3 , pm(0, 1,1,1) = ¥>m(1,0,0,0) = & 4 , 
It is shown in Figure [21] how vertices of the unit tesseract in Z are mapped by 




Figure 22. Mapping 92^ on the unit tesseract. 



It is straightforward to prove the following 
Proposition 4.4. ipM is a double reflection net. 

Poncelet-Darboux grids and double reflection nets. Let £ be an ellipse in the Eu- 
clidean plane: 

2 2 
x v 
8 : — + a > 6 > 0, 

a 

and (ai)igz a billiard trajectory within £. 

As it is well known, all lines <Zj are touching the same conic C confocal with £ . 
Here, we will additionally suppose that C is an ellipse. Denote by A the set of 
tangents of C. 

Fix m positive integers k\, . . . , k m and define the mapping: 

f D : Z m ->• A, ip D (ni, . . . ,n m ) = a„ lkl+ ... +nmkm . 
Proposition 4.5. Map tpo is a double reflection net. 
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Proof. Si nce (fD{n + iej) = a ni k 1 +-+n m k m +ik :i , (n = (n%, . . . , n m )), it follows by 
|DR201l| Theorem 18] that sequence (<^£>(n + ie.j))ie7. represents a billiard trajec- 
tory within some ellipse £j, confocal with £ and C. 

Immediately, by Definition ^. 11[ lines <f(n ), (/?(n +ej), ip(n +ej), ip(n +ei+ej) 
form a virtual reflection configuration for each no G Z m , i, j G {1, . . . , m}. 

Moreover, by Propositon 12.131 they also form a double reflection configuration. 

□ 

Remark 4.6. It is interesting to consider only nets where m ellipses £j appearing 
in the proof of Proposition \4-5\ are distinct. If some of them coincide, then we may 
consider a corresponding subnet. 

Suppose that (a,) is a non-periodic trajectory. Then, choosing any m, and any 
set of distinct positive numbers kt, k m , we get substantially different double 
reflection nets. 

For {fli) being n-periodic, it is enough to consider the case fej = i, i S {1, . . . , [n/2]}, 
(m = [n/2]). 

Example 4.7. Suppose (a,) is a 5-perodic billiard trajectory within £, see Figure 




Figure 23. A Poncelet pentagon. 

The corresponding double reflection net is: 

ip D : Z 2 -> A, y£)(ni,n 2 ) = a ni+2 „ 2 . 

s-skew lines and double reflection nets. Now, let us consider a family of confocal 
quadrics in E d (d > 3) and fix its d — 1 quadrics. As usually, A is the set of all 
lines tangent to the fixed quadrics. 

It is shown in [DR2008 that, from a line in A, we can obtain any other line from 
that set in at most d—1 reflections on quadrics from the confocal family. We called 
lines a, b from A s-skew if s is the smallest number such that they can be obtained 
by s + 1 such reflections. 

Now, suppose lines a, b are s-skew (s > 1), and let Qi, . . . , Qs+i be the corre- 
sponding quadrics from the confocal family. 
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Theorem 4.8. There is a unique double reflection net 

ip s : Z S+1 ^A 



which satisfies the following: 

• <p s (0, ... ,0) = a; 

• <£>s(l) ••..!) = &; 

• {tp(no + zej)}jgz represents a billiard trajectory within Qj, for each j G 
{1, . . . ,s + 1} and n G Z s+1 . 

Proof. First, using the construction of the billiard algebra from [DR2 008] . we are 
going to define mapping ip s on {0, 1} S+1 . 

For a permutation p = (pi, . . . ,p s +i) of the set {f , . . . , s + 1}, we take a sequence 
of lines (^q, . . . ,£f +1 ) such that £q = a, £^ — b, and if satisfy the reflection 
law off Q Pi for each i G {f , . . . , s + 1}. Such a sequence exists and it is unique. 
Moreover, if k G {f, . . . ,s + 1} is given, and permutations p, p' coincide in the 
first k coordinates, then If = if for i < k. Take {ii,...,ifc} to be a subset 
of {1, . . . , s + 1 }, and p any permutation of set {1, . . . , s + 1} with pi = i\, . . . , 
Pk = ik- We define: 



where x — X{ii, ...,-**,} is the corresponding characteristic function on {f , . . . , s + 1}: 



In this way, we constructed ip s on {0, 1} S+1 . 

Subsequently, <p s can be extended to the rest of Z s+1 , so that {</? s (no + iej)},- e z 
will represent billiard trajectories within Qj. 



Construction of double reflection nets. Let Qi, Q m be distinct quadrics be- 
longing to a confocal family and i a line in P d . Let us choose lines £, satisfying 
with I the reflection law off Q,, 1 < i < m. 

Theorem 4.9. There is a unique double reflection net Lp : Z" 1 — > At, with the 
following properties: 



• ip(ei) — ii, for each i G {I, . . . ,m}. 

By Ai, we denoted the set of all lines in P d touching the same d — 1 quadrics from 
the confocal family as i. 

Proof. First, we define tp on {0, 1}"\ from the condition that lines corresponding 
to each 2-face of the unit cube need to form a double reflection configuration. 
This construction is unique because of Proposition 12.141 and correct, due to the 
3£>-consistency property proved in Theorem 14.21 

At all other points of Z m , ip is uniquely defined from the request that {<^(no + 
iej)}i£z will be billiard trajectories within Qj. 

Consistency of the construction follows again from Thcorcm l4.2l □ 




This construction is correct and unique due to Theorem 14.21 



□ 
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Focal nets and F -transformations of double reflection nets. Let tp : Z m — > A be a 
double reflection net. 

For given n G Z m and distinct indices k G {1, . . . m}, consider the following 
points of its i-th focal net: 

Fi = F (l \n a ) = </j(n ) n ip(n + e t ), 

F i3 = F {l) {n Q + e,-) = tp(n + e,-) n <^(n + e 3 - + e,), 

F ik = F^(n + efc) = ^(n + e fe ) n <ys(n + e fe + e.,), 

F ijk = F (l) (n + ej + efe) = <^(n + ej + e k ) n ^(n + ej + e k + e»). 

Proposition 4.10. Points Fi, Fy, Fik, -Pijfc are coplanar. 

Proof. This is a consequence of the theorem of focal nets from BS2008 . However, 
we will show the direct proof, from configurations considered in Section l4.ll 

The four points belong to quadric Q4. The tangent planes to Q4 at these points, 
divided into two pairs, determine two pencils of planes. According to Theorem 
14.11 the two pencils are coplanar; thus they intersect. As a consequence, the lines 
of poles with respect to the quadric Qi, which correspond to these two pencils of 
planes, also intersect. It follows that the four points are coplanar. □ 

We are going to construct an ^-transformation of the double reflection net. 

First, we select a quadric Q$ from the confocal family and introduce line £' which 
satisfies with <^(no) the reflection law on Qg. 

By Theorem l4.9[ it is possible to construct a double reflection net (p : Z m+1 — > 
A, such that: 

• £(n,0) = <p(n); 
. £(n ,l)=l'. 

Now, we define: 

ip+ : Z m ^A e , p+(n)=$j(n,l). 
Proposition 4.11. Map ip + is an F -transformation of tp. 

Proof. Lines y + (n) and <p(n) intersect, since they satisfy the reflection law off 
Qs- □ 

Double reflection nets and Grassmannian Darboux nets. Let us recall the definition 
of a Grassmannian Darboux net from [ABS2009 : a map from the edges of a regular 
square lattice Z m to the Grassmannian Gf. of r-dimensional projective subspaces 
of the d-dimensional projective space is a Grassmanian Darboux net if the four 
r-spaces of an elementary quadrilateral belong to a (2r + l)-space. For r = 0, 
the ordinary Darboux nets from |Sch2003] are obtained, where the four points of 
intersection associated to a quadrilateral, belong to a line. 
Now, consider a general double reflection net (|4.ip . 

To each edge (no, no + e,-) of Z m , we can associate the plane which is tangent to 
Qi at point <^(no) n <p(rio + ej). Since the lines corresponding to the vertices of a 
face form a double reflection configuration, the four planes associated to the edges 
belong to a pencil. 
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In this way, we see that a double reflection net induces a map: 

E(Z m ) -> GU, 

where E(Z m ) is the set of all edges of the integer lattice Z m . 

In this way, double reflection nets induce a subclass of dual Darboux nets. 

It was shown in |Sch2003] how to associate discrete integrable hierarchies to the 
Darboux nets. 

4.3. Yang-Baxter map. A Yang-Baxter map is a map R : X x X X x X, 

satisfying the Yang-Baxter equation: 

-R23 R13 R12 — R12 R13 R23, 

where Rij : X X X X X — > X X X X X acts as R on the i-th and j-th factor in the 
product, and as the identity on the remaining one, see [ABS2004] and references 
therein. 

Here, we are going to construct an example of Yang-Baxter map associated to 
confocal families of quadrics. To begin, we fix a family of confocal quadrics in CP": 

< 4 - 2 > * = 54j t + - + s^r-^ 

where ai, . . . , ad are constants in C, and [z% : Z2 '■ • ■ ■ '■ z n +i] are homogeneous 
coordinates in CP". 

Take X to be the space CP"* dual to the n-dimensional projective space, i.e. the 
variety of all hyper-planes in CP" . Note that a general hyper-plane in the space 
is tangent to exactly one quadric from family (|4.2|) . Besides, in a general pencil of 
hyper-planes, there are exactly two of them tangent to a fixed general quadric. 

Now, consider a pair x, y of hyper-planes. They are touching respectively unique 
quadrics Q a , Qp from (|4.2I) . Besides, these two hyper-planes determine a pencil of 
hyper-planes. This pencil contains unique hyper-planes x', y' , other than x, y, that 
are tangent to Q a , Qp respectively. 

We define R : CP"* x CP"* CP"* x CP"*, in such a way that R(x,y) = 
(x',y') if (x',y') are obtained from (x,y) as just described. 

Maps 

R12, Ria, R23 ■ CP"* x CP"* x CP"* -> CP"* x CP"* x CP"* 
are then defined as follows: 

Ri2(x,y,z) = {x',y',z) for (V, y') = R(x, y); 

R-Ls(x,y,z) = (x',y,z') for (x',z') = R(x, z); 

R-zs{x,y,z) = (x,y',z') for (y' , z') = R(y , z) . 

To prove the Yang-Baxter equation for map R, we will need the following 

Lemma 4.12. Let Q a , Qp, Q 7 be three non- degenerate quadrics from family \4-<ty 
and x, y, z respectively their tangent hyper-planes. Take: 

(x2,yi) = R{x,y), (x s ,zi) = R(x,z), (y 3 ,z 2 ) = R(y,z). 

Let X23, J/13, 212 be the joint hyper-planes of pencils determined by pairs (2:3, 2/3) 
and {x2,Z2), (£3,2/3) and (yi,zi), (yi,zi) and (£2,^2) respectively. 
Then X23, J/13, Z12 touch quadrics Q a , Qp, Q 7 respectively. 
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Proof. This statement, formulated for the dual space in dimension n = 2 is proved 
as [ABS2004( Theorem 5]. 

Consider the dual situation in an arbitrary dimension n. The dual quadrics Q* , 
Q%, Qly belong to a linear pencil, and points x* , y* , z* , dual to hyper-planes x, 
y, z, are respectively placed on these quadrics. Take the two-dimensional plane 
containing these three points. The intersection of the pencil of quadrics with this, 
and any other plane as well, represents a pencil of conies. Thus, Theorem 5 from 
ABS2004J will remain true in any dimension. 

This lemma is dual to this statement, thus the proof is complete. □ 

Theorem 4.13. Map R satisfies the Yang-Baxter equation. 

Proof. Let x, y, z be hyper-planes in CP™. We want to prove that 

R23 o R13 o Ru(%, y, z) = i?i2 o R13 o Rn{x, y, z). 

Denote by Q a , Qp, Q 1 the quadrics from (|4. 2[) touching x, y, z respectively. 
Let: 

(X,y,Z) ► [X 2 ,yi,Z) > (X23,yi,Zi) > (X23,yi3,Zi2), 

(x,y,z) > (x,y 3 ,z 2 ) > {x 3 ,y 3 ,z 12 ) > {x 23 , y 13 , z 12 ). 

Now, apply Lemma 14.121 to hyper-planes x, y, z 2 . Since: 

(x 2 ,yi) = R(x,y), (x 3 ,z' 12 ) = R(x,z 2 ), (y 3 , z) = R(y, z 2 ), 

we have that the joint hyper-plane of pencils (x 3 , y 3 ) and (x 2 , z) is touching Q a - 
therefore, this plane must coincide with x 23 and x' 23 , i.e. 0:23 = x' 23 . Also, the joint 
hyper-plane of pencils (y\, z' 12 ) and (x 2 , z) is touching Q 7 - therefore, this is z\ and 
Z\ 2 = z' 12 . Finally, the joint hyper-plane of pencils (x 3 , y 3 ) and (yi,z' l2 ) is tangent 
to Q/3 - it follows this is yi 3 — y' 13 , which completes the proof. □ 

Remark 4.14. Instead of defining R to act on the whole space CP n * x CP™*, we 

can restrict it to the product of two non-degenerate quadrics from j4.%fy , namely: 

R(a,[3) : Q* x Q£ -> Q* x Q% 

where pair (x,y) of tangent hyper-planes is mapped into pair (x\,y\) in such a way 
that x, y, x\, y\ belong to the same pencil. 
The corresponding Yang-Baxter equation is: 

#23(^,7) #13 («, 7) Ri2(a,P) = Ri 2 {a,P) o R 13 (a,<y) o R 23 (a,f3), 
where both sides of the equation represent maps from Q* a x Q*j x Q* to itself. 

In [ABS2004J, for irreducible algebraic varieties X\ and X 2 , a quadrirational 
mapping F : X\ x X 2 is defined. For such a map F and any fixed pair (x,y) € 
X\ x X 2 , except from some closed subvarieties of codimension at least 1, the graph 
Pf c X\ x X 2 x X\ x X 2 intersects each of the sets {x} x {y} x X\ x X 2l X\ x X 2 x 
{x} x {y}, X\ x {y} x {x} x X 2 , {x} x X 2 x X\ x {y} exactly at one point (see 
ABS2004. Definition 3]). In other words, Tf is the graph of four rational maps: 
F, F- 1 , F, F- 1 . 

The following Proposition is a generalization of [ABS2004, Proposition 4]. 
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Proposition 4.15. Map R(a,(3) : Q* x Q*^ -» Q* x Ql, is quadrirational. It is 
an involution and it concides with its companion R(ct,j3). 

5. PSEUDO-INTEGRABLE BILLIARDS AND LOCAL PONCELET THEOREM 

5.1. Billiards in domains bounded by a few confocal quadrics. In this sec- 
tion we analyze billiard dynamics in a domain bounded by arcs of a few confocal 
conies, such that there are reflex angles on the boundary. In order to describe some 
phenomena appearing in such systems, let us consider the domain T>q bounded by 
two confocal ellipses from family (|2.ip and two segments placed on the smaller axis 
of theirs, as shown in Figure [24j More precisely, we fix parameters fa, fa such that 




Figure 24. Domain bounded by two confocal ellipses and two 
segments on the y-axis. 

fa < fa < b, and take the border of T>o to be: 

dv = r 1 ur 2 ur 3 ur 4 , 

T 1 ={(x,y)eC l3l \x>0}, 

r 2 = {(x,y)eCp 2 \x<0}, 
r 3 = {(0,y) \Vb^P~2<y<Vb~Ti}, 

r 4 = {(o, y ) 1 -Vb-Pi<y<-Vb-fo}- 

Notice that segments T^, T4 are lying on the the degenerate conic C a of family (|2.1I) . 

By Chasles'theorem Chal827 , as we have already mentioned in Section 12. 1[ 
each line in the plane is touching exactly one conic from a given confocal family 
- moreover, this conic remains the same after the reflection on any conic from 
the family. Thus, each billiard trajectory in a domain bounded by arcs of several 
confocal conies has a caustic from the confocal family. 

Consider billiard trajectories within domain T>o whose caustic is an ellipse C ao 
completely placed inside the billiard table, i.e. fa < «o < b. An example of such a 
trajectory is shown in Figure [231 

Such billiard trajectories fill out the ring 1Z placed between the billiard border 
and the caustic, see Figure |2"o1 
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Figure 25. A billiard trajectory in T>q with an ellipse as caustic. 




Figure 26. Ring K. 

Let us examine the leaf of the phase space composed by these trajectories. This 
leaf is naturally decomposed into four rings equal to 1Z, which are glued with each 
other along the border segments. Let us describe this in detail: 

IZi This ring contains the points in the phase space that correspond the billiard 
particle moving away from the caustic and the clockwise direction around 
the ellipses center. 

IZ2 Corresponds to the motion away from the caustic in the counterclockwise 
direction. 

IZ3 Corresponds to the motion towards the caustic in the counterclockwise 
direction. 

1Z± Corresponds to the motion towards the caustic in the clockwise direction. 

Let us notice that the reflection off the two ellipse arcs contained in the billiard 
boundary changes the direction of the particle motion with respect to the caustic, 
but preserves the direction of the motion around the foci. The same holds for 
passing though tangency points with the caustic. On the other hand, reflection 
on the axis changes the direction of motion around the foci, but preserves the 
direction with respect to the caustic. Thus, the four rings are connected to each 
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other according to the following scheme: 



n 2 




Let us represent all the rings in Figures [57] and 




T4 T4 1^4 r 4 

Figure 27. Rings Hi, K 2 , K 3 , K 4 . 
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Figure 28. Gluing rings TZ U TZ 2 , K 3 , Tl A . 



Now, we have the following 

Proposition 5.1. All billiard trajectories within domain T>o with a fixed elliptical 
caustic form an orientable surface of genus 3. 

5.2. Examples. In this section, we are going to analyze examples of billiards in 
domains described in Section 15.11 in a special case when confocal family (|2.1[) is 
degenerate: a = b. The family then consists of concentric circles. It is convenient 
to analyze such cases, since they may be approached with elementary means, yet 
all phenomena appearing in non-degenerate cases are to be preserved, as will be 
shown in Section 15.51 
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Domain bounded by circles with rotation numbers — and — . Let us consider example 
y 3 4 v 

of the billiard within the domain with two concentric half-circles of radii 2R, Ry/2 
and the corresponding segments. We will consider caustic with radius R. 

In this case, there exist six trajectories connecting singular points corresponding 
to reflex angles of the boundary. Those trajectories are represented in Figure 1291 
Each polygonal line shown on the figure corresponds to two trajectories in the phase 
space, depending on direction of the motion. 




Figure 29. Saddle-connections corresponding to circles with ro- 
tation numbers — and — . 

3 4 

Vertices of the saddle-connections divide the billiard border into thirteen parts, 
see Figure[30l For the definition of saddle-connections - see Section[T] Introduction. 



Figure 30. Parts of the boundary corresponding to circles with 

rotation numbers — and — . 

3 4 

All trajectories in this billiard domain corresponding to the fixed caustic are 
periodic: 

• either all bouncing points of a given trajectory are in gray parts - in this 
case the billiard particle hits twice each gray part until the trajectory be- 
comes closed. Such a trajectory is 12-periodic, with with four bounces on 
the smaller circle, six bounces on the bigger one, and one on each of the 
segments on the y-axis (see Figures l3lk and l3"Ib ) ; 
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• or all bouncing points are in orange parts - the particle will hit each part 
once until closure and the trajectory is 7-periodic. Such a trajectory hits 
three times the bigger circle and four times the smaller one (see Figure 

[Mb). 




The corresponding level set in the phase space is divided by the saddle-connections 
into three parts: 

• the part containing all 12-periodic trajectories: this part is bounded by four 
saddle-connections whose projections on the configuration space is shown 
in Figures l29b and [29b; 

• two parts containing all 7-periodic trajectories winding about the caustic in 
the clockwise and counterclockwise direction: these parts are bounded by 
saddle-connections winding in the same direction whose projections on the 
configuration space is shown in Figures l2"9"b andl2"9"b; the saddle-connections 
corresponding to Figure l2"9"b are lying within the corresponding parts. 

Domain bounded by circles with rotation numbers — and — . Let us consider example 
... . ..46 . . 

of the billiard within the domain determined with two concentric half-circles with 

rotation numbers equal to — and — . 

4 6 

In this case, there exist six saddle-connections, represented in Figure [321 Each 
polygonal line shown on the figure corresponds to two trajectories in the phase 
space, depending on direction of the motion. 

Vertices of the saddle-connections divide the billiard border into eight parts, see 
Figure [Ml 

All trajectories in this billiard domain corresponding to the fixed caustic are 
periodic: 

• either all bouncing points of a given trajectory are in gray parts - in this case 
the billiard particle hits twice each gray part until the trajectory becomes 
closed and the trajectory is 6-periodic, see Figures 1Mb and [Mb: 

• or all bouncing points are in orange parts - the particle will hit each part 
once until closure and the trajectory is 5-periodic, see Figure [Mb. 
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Figure 32. Saddle-connections corresponding to circles with ro- 
tation numbers — and — . 

4 6 




rotation numbers — and — . 

4 6 




Figure 34. Periodic trajectories corresponding to circles with ro- 
tation numbers — and — . 

4 6 



The corresponding level set in the phase space is divided by the saddle-connections 
into three parts: 
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• the part containing all 6-periodic trajectories: this part is bounded by four 
saddle-connections whose projections on the configuration space is shown 
on Figures 152k and 152b; 

• two parts containing all 5-periodic trajectories winding about the caustic in 
the clockwise and counterclockwise direction: these parts are bounded by 
saddle-connections winding in the same direction whose projections on the 
configuration space is shown in Figures 132k and 132b; the saddle-connections 
corresponding to Figure [3"2"b are lying within the corresponding parts. 

5 — \/5 y/5 

Domain bounded by circles with rotation numbers and — . In this example, 

. . 10 10 

there exist six saddle-connections, represented in Figure 1351 Each polygonal line 

shown on the figure corresponds to two trajectories in the phase space, depending 

on direction of the motion. 




Figure 35. Saddle-connections corresponding to circles with ro- 

i 5-V5 y/E 

tation numbers and — . 

10 10 



Vertices of the saddle-connections divide the billiard border into eleven parts, 
see Figure [36] 



Figure 36. Parts of the boundary corresponding to circles with 

rotation numbers ^— and ^— . 

10 10 

All trajectories in this billiard domain corresponding to the fixed caustic are 
periodic: 
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• either all bouncing points of a given trajectory are in gray parts - in this case 
the billiard particle hits twice each gray part until the trajectory becomes 
closed and the trajectory is 14-periodic, see Figures 157k and [57b . Notice 
that such a trajectory bounces six times on each of the circles and once on 
each of the segments; 

• or all bouncing points are in orange parts - the particle will hit each part 
once until closure and the trajectory is 4-periodic, see Figure [37b. Such a 
trajectory reflects twice on each of the circular arcs. 




Figure 37. Periodic trajectories corresponding to circles with ro- 
tation numbers — anc ^ "Jq"' 

The corresponding level set in the phase space is divided by the saddle-connections 
into three parts: 

• the part containing all 14-periodic trajectories: this part is bounded by four 
saddle-connections whose projections on the configuration space is shown 
on Figures [35k and [35k . The saddle-connections corresponding to Figure 
[35b are lying inside this part; 

• two parts containing all 4-periodic trajectories winding about the caustic 
in the clockwise and counterclockwise direction: these parts are bounded 
by saddle-connections winding in the same direction whose projections on 
the configuration space is shown in Figures [35k and 135b . 

5.3. Topological estimates. Let T> be a bounded domain in the plane such that 
its boundary T = dT> is the union of finitely many arcs of confocal conies from the 
family 

We consider the billiard system within T>. Any trajectory of this billiard will 
have a caustic - a conic from (|2.ip touching all lines containing segments of the 
trajectory. Let us fix C\ as caustic. 

Notice that all tangent lines of a conic fill out an infinite domain in the plane: if 
the conic is an ellipse, the domain is its exterior; for a hyperbola, it is the part of 
the plane between its branches. 

Denote by V\ the intersection of T> with the domain containing tangent lines 
of caustic C\ . All billiard trajectories with caustic C\ are placed in V\ Q . T>\ is 
a bounded set whose boundary T\ = dV\ is the union of finitely many arcs of 
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conies from (|2.1[) . We assume that V\ is connected as well, otherwise we consider 
its connected component. 

All billiard trajectories in domain T> with the caustic C\ will correspond to a 
certain compact leaf M.(\ ) in the phase space. M.\ is obtained by gluing four 
copies of V\ along the corresponding arcs of the boundary Y\ a = dT>\ , similarly 
as it is explained in Section [Ol 

On M.\ , singular points of the billiard flow correspond to vertices of reflex 
angles on the boundary of T>\ . Since confocal conies are orthogonal to each other 
at points of intersection, only two types of such angles may appear: full angles and 
angles of 270°. A vertex of a full angle is the projection of two singular points in the 
phase space. Each of the singular points has four separatrices. On the other hand, 
a vertex of a 270° is a projection of only one singular point having six separatrices. 

Using the Euler-Poincare formula, as in [Via2008 , we get the following estimate 
for the total number N — N(J\4\ ) of saddle-connections: 

Proposition 5.2. The total number N — N(A4\ ) of saddle- connections is bounded 
from above: 

k 

N(M Xo )<~J2 s i = k-x{Mx ), 
i=i 

where k is the number of singular points of the flow on M\ , and s\, Sk 
numbers of separatrices at each singular point. 

As a corollary, we get the following 

Proposition 5.3. Consider billiard within T> with C\ as a caustic. If the corre- 
sponding subdomain T>\ has k reflex angles on its boundary T\ then: 

• N < 3k; 

• g(M Xo ) = k + l. 

Notice that the genus of the surface M.\ depends only on the number of reflex 
angles on the boundary of T>\ and not of their types. Also, k < k. 

Example 5.4. • If there are no reflex angles on the boundary, i.e. k = 0, 

then M\ is a torus: g = 1, N = 0; 

• if there is only one reflex angle on the boundary, independently if it is a 
270° angle or a full angle, we have that g = 2. 

It is well known that the Liouville- Arnold theorem (see |Arnl978] 'l describes 
regular compact leaves of a completely integrable Hamiltonian system as tori, with 
dymanics being quasi-periodic on these inavriant tori. On some of the tori, the 
dynamics is exclusively periodic and for a fixed such a torus, the period is fixed. 

We finish this section by formulating of an analogue of the Liuoville- Arnold 
theorem for pseudo-integrable billiard systems. It is a consequence of the Maier 
theorem from the theory of measured foliations (see [Mail943 ( Via2008 ). In our 
case, the measured foliation is defined by the kernel of an exact one- form (3\ := dK\, 
where functions K\ are defined in Section |2~T1 

Theorem 5.5. There exist paiwise disjoint open domains D\, D n on A^Aq, 
each of them being invariant under the billiard flow, such that their closures cover 
M.\ a and for each j € {1, . . . , n}: 
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• either Dj consists of periodic billiard trajectories and is homeomorphic to 
a cylinder; 

• or Dj consists of non-periodic trajectories all of which are dense in Dj . 
The boundary of each Dj consists of saddle-connections. 

We see, that in contrast to completely integrable Hamiltonian systems, compact 
leaves of our billiards could be of a genus greater than 1. Moreover, one leaf could 
contain regions with periodic trajectories with different periods for different regions 
and simultaneousely could contain regions with non-periodic motion. Because of 
that, we call such systems pseudo-integrable, taking into account the fact that they 
possess two independent commuting first integrals, as it has been shown in Section 

5.4. Poncelet theorem and Cayley-type conditions. For billiards within con- 
focal conies without reflex angles on the boundary, it is well known that the famous 
Poncelet porism holds (see |DR2004ilDR2011j l: 

(A) if there is a periodic billiard trajectory with one initial point of the boun- 
dary, then there are infinitely many such periodic trajectories with the same 
period, sharing the same caustic; 

(B) even more is true, if there is one periodic trajectory, then all trajectories 
sharing the same caustic are periodic with the same period. 

However, when reflex angles exist, which is the case studied in the present paper, 
one can say that (A) is still generally true. However, (B) is not true any more. In 
other words, the Poncelet porism is true locally, but not globally. 

Theorem 5.6. There exist subsets 5\, . . . , S n of the boundary T\ , with the fol- 
lowing properties: 

• Si, . . . , S n are invariant under the billiard map; 

• 8\, . . . , 5 n are pairwise disjoint; 

• each Si is a finite union of di open subarcs of T\ : 

[J r r 

3=1 

• closure of Si U • • • U Sn is T, 
such that they satisfy: 

• if one billiard trajectory with bouncing points within Si is periodic, then all 
such trajectories are periodic with the same period m. Moreover, m is a 

Ui 

multiple of di and every such a trajectory bounces the same number — of 

d l 

times off each arc ; 

• if billiard trajectories having vertices in Si are non-periodic, then the bounc- 
ing points of each trajectory are dense in Si. 

The boundary of each Si consists of bouncing points of saddle-connections. 

This theorem is a consequence of Theorem 15.51 from the previous section. The 
proof follows from the fact that each of the domains Di intersects the boundary 
T\ and forms Si — T\ n Di. 
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In |DR2004j condit ions of Cayley's type for periodicity of billiards within several 
confocal quadrics in the Euclidean space of an arbitrary dimension were derived, 
see also DR2006a for detailed examples. 

We analyzed there billiards within domains bounded by arcs of several confocal 
quadrics and the billiad ordered game within a few confocal ellipsoids. Unlike in 
the present article, domains considered in [DR2004 ( DR2006a did not contain reflex 
angles at the boundary. However, the technique used there to describe periodic 
trajectories can be directly transferred to the present problems. 

Before stating the Cayley-type conditions, recall that a point is being reflected 
off conic C\ from outside if the corresponding Jacobi elliptic coordinate achieves 
a local maximum at the reflection point, and from inside if there the coordinate 
achieves a local minimum (see [DR2004] ). 

Theorem 5.7. Consider domain T> bounded by arcs of k ellipses Cp 17 . . . , Cp k , I 
hyperbolas C 7l , . . . , C 7i , and several segments belonging to degenerate conies from 
the confocal family \2.1\) : 

Pi,..., 0k G (-00, b), k> 1, 71,..., 7; G (6, a), I > 0. 

Let C ao be an ellipse contained within all ellipses Cp 1 , . . . , Cp k : b > ao > Pi for 
all i G {1, . . . , k}. A necessary condition for the existence of a billiard trajectory 
within T> with C ao as a caustic which becomes closed after: 

• n[ reflections from inside and n'[ reflections from outside offC^, 1 < i < k; 

• m'j reflections from inside and m'j reflections from outside off C 7 . , 1 < j < 
I; 

• total number of p intersections with the x-axis and reflections off the seg- 
ments contained in the x-axis; 

• total number of q intersections with the y-axis and reflections off the seg- 
ments contained in the y-axis; 

is: 

k I 

J2« - <)M(P ft ) - A(P ao )) + £(m$ - m'!)A(P^) + P A(P a ) - qA(P b ) = 0, 
i=i 3=1 

m'j — m'j + p — q = 0. 
Here A is the Abel-Jacobi map of the ellitic curve: 

T : s 2 =V(t) := (a-t)(b-t)(a -t), 
and P$ denotes point (5, yJV{5)) on T . 

Proof. Following Jacobi [Jacl884] and Darboux |Darl870| . similarly as in |DR2004] . 
we consider sums 

r dXi r d\ 2 and r XidM f A 2 rfA 2 

over billiard trajectory A\ . . . An ■ Here (Ai,A2) are Jacobi elliptic coordinates, 
Ai < A 2 . The second integral is equal to the length of the trajectory, while the first 
one is zero. 

Notice that, along a trajectory, Ai achieves local extrema at points of reflection 
off ellipses and touching points with the caustic, and A2 at points of reflection 
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off hyperbolas and intersection points with the coordinate axes, we obtain that 
A\ = An is equivalent to the condition stated. □ 

We illustrate this theorem on the example when the billiard table is T>o, as 
defined in Section |5~T1 

Example 5.8. A necessary condition for the existence of a billiard trajectory within 
T>o withC ao as a caustic, such that it becomes closed after n\ reflections offCp 1 and 
712 reflections offCp 2 is: 

niA(P Pl ) + n 2 A{P h ) = {n 1 +n 2 )A{P ao ). 

Notice that in this case number p and q are always even and equal to each other. 
Since 2A(P a ) = 2A(Pb), the corresponding summands are cancelled out. 

5.5. Interval exchange transformation. In this section, we are going to estab- 
lish a connection of the billiard dynamics within domain 2?o defined in Section 15.11 
with the inteval exchange transformation. 

Interval exchange maps. Let I C R be an interval, and {I a \ a € A} its finite 
partition into subintervals. Here A is a finite set of at least two elements. We 
consider all intervals to be closed on the left and open on the right. 

An interval exchange map is a bijection of / into itself, such that its restriction 
on each I a is a translation. Such a map / is determined by the following data: 

• a pair {ttq,ttx) of bijections A — > {1, ...,d} describing the order of the 
subintervals {I a } in / and {/(/<*)} m /(-0 = !• We denote: 

/ 7^(1) 7T - 1 (2) ... ^\d) \ 

v Trf^l) 7r 1 - 1 (2) ... n^(d) )■ 

• a vector A = (A q ) q6 ^i of the lengths of I a . 

Billiard dynamics. To each billiard trajectory, we join the sequence: 

where X n are joint points of the trajectory with the caustic, while s n — + if at X n 
the trajectory is winding counterclockwise and s„ = — if it is winding clockwise 
about the caustic. 

Introduce metric \x on the caustic C ao as in Proposition 12.11 Then, we parame- 
trize C ao by parameters: 

p : C ao -> [0,1), q : C ao -+ [-1,0), 

which arc natural with respect to /i such that p is oriented counterclockwise and q 
clockwise along C ao , and the values p = and q = — 1 correspond to points Pq, Qo 
respectively, as shown in Figure 

Consider one segment of a billiard trajectory, and let X G C ao be its touching 
point with the caustic. Suppose that the particle is moving counterclockwise on 
that segment. From Figure [551 we conclude: 

• if X is between points Pi and Pi then the particle is going to hit the arc 

• if X is between P2 and Pq, the particle is going to hit the arc C\ 1 ; 
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• for X between Pq and Pi, the particle is going to hit and the upper 
segment before the next contact with the caustic and the direction of motion 
is changed to clockwise. 

Similarly, if the particle is moving in clockwise direction, we have: 

• if X is between points Qi and Q 2 then the particle is going to hit the arc 

C\ 2 ] 

• if X is between Q 2 and Qq, the particle is going to hit the arc C\ 1 ; 

• for X between Q and Qi, the particle is going to hit C\ 1 and the lower 
segment before the next contact with the caustic and the direction of motion 
is changed to counterclockwise. 

To see the billiard dynamics as an inteval exchange transformation, we make the 
following identification: 

(*,+)~p(X), (X,-)~ q (X). 

In other words: 

• we identify the joint point X of a given trajectory with the caustic with 
p(X) € [0, 1) if the particle is moving in the counterclockwise direction on 
the corresponding segment; 

• for the motion in the clockwise direction, we identify X with q(X) S [—1,0). 

Denote the rotation numbers ri = p(Xi), r 2 = p(A 2 ) (see Proposition ^. ip . 
The parametrizations values for points denoted in Figure [38] are: 

p(P o ) = 0, p(Pi) = r 1 -r 2 , p(P 2 )=ri-r 2 + ~, 

q{Qo) = -1, q(Qi) = n - r 2 - l, qiQi) = n -r 2 - -. 

Now, we distinguish three cases depending on the position of point P with 

1 r 2 

respect to the x-axis (see Figure [35]) , i.e. on the sign of - + — — r x . 
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Pq is on the x-axis: — I n =0. The interval exchange map is: 
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[i-n,i), 



as shown in Figure 1391 

AS C D F F 

C B D F E A 

1 ^\\\\\\\\\\\vf^66666666q , 



Figure 39. Interval exchange transformation for the case 4 + ^ — 

n = 0. 



To the map, pair (7T, A) is joined: 

A B C D E F 
^ C B D F E A 

if 1 11 1 

1 r 2 

Pq is above the x-axis: — + — — rj > 0. The interval exchange map in this case is 
shown in Figure l40l and given by: 
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Figure 40. Interval exchange transformation for the case \ + ^ — 

n > 0. 
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The map can be desribed by the pair (tt, A): 

A B C D E F G 
D B E C H F A 



H 

G 



n 



2ri + -, r l5 ri 



1 



T2, ^, r 2 



2n 



2' ri 



Po *s below the x-axis: 
to the billiard dynamics is: 
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ri < 0. The interval exchange map corresponding 
-l,-|-rO 

-5 - n>n. - ^2 - 1) 

Ti - r 2 - l,ri - r 2 - 5) 
h-ra-1,0) 

h-ri) 
- n,n - r 2 ) 
'"1 -r 2 ,ri -r 2 + |) 



ri - r 2 + n> l)j 




Figure 41. Interval exchange transformation for the case x + ^r 

n < 0. 

To the map, pair (7r, A) is joined: 

(\ 1111 III 

A=(j-n, 2rx-r 2 --, 3,^ + --^ --n, 2n - r 2 - -, ^r 2 + --n 

Notice that in all three cases the interval exchange transformations depend only 
on the rotation numbers r%, r 2 . Thus, we got 

Theorem 5.9. The billiard dynamics inside the domain T>o with ellipse C ao as the 
caustic, does not depend on the parameters a, b of the confocal family but only on 
the rotation numbers r\, r 2 . 

5.6. Keane condition and minimality. An interval exchange transformation is 
called minimal if every orbit is dense in the whole domain. When considering 
pseudo-billiards, minimal interval exchange transformations will correspond to the 
cases when all orbits are dense in the domain between the billiard border and the 
caustic. 
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Following Via2008, Zor2006], we are going to formulate a sufficient condition for 
minimality. Let / be an interval exchange transformation of /, given by pair (jr, A). 
Denote by p a the left endpoint of I a . Then the transformation satisfies the Keane 
condition if: 

f m (p a ) + Pp for aU m > 1, a G A, P G A \ {^(l)}. 



Obviously, none of the transformations from Section 15.51 satisfies the Keane condi- 
tion: namely, the midpoint of the interval is the left endpoint of one of I ai and it 
is the image of another endpoint in the corresponding interval exchange map. 

The goal of this section is to find an analogue of the Keane condition for interval 
exhange transformations appearing in the billiard dynamics. 

Billiard-like transformations and modified Keane condition. Analysis of the exam- 
ples from Section l5~5l motivates the following definitions. 

Definition 5.10. An interval exchange transformation f of I = [—1, 1) is billiard- 
like if the partition into subintervals satisfies the following: 

• for each a, I a is contained either in [—1,0) or [0, 1); 

• for each a, f(I a ) is contained either in [—1,0) or [0, 1); 

• both [—1,0) and [0, 1) contain at least two intervals of the partition. 

Definition 5.11. We will say that a billiard-like interval exchange transformation 
f satisfies the modified Keane condition if 

f m (pu) 7^ P/3 for all m>\, a G A, and (3 G B such that pp ^ {—1,0}. 

Lemma 5.12. If a billiard-like interval exchange transformation satisfies the mod- 
ified Keane condition, then the transformation has no periodic points. 

We say that an interval exchange transformation is irreducible if for no k < \A\ 
the union 

T Q (1) T Q O) 

is invariant under the transformation. The usual Keane condition implies irre- 
ducibility. However, this is not the case for the modified Keane condition - it may 
happen that the transformation falls apart into two irreducible transformations on 
[— 1, 0) and [0, 1). On the other hand, if for a transformation satisfying the modi- 
fied Keane condition there is an interval I a C [—1,0) such that f(I a ) C [0, 1), the 
irreducibility will also take place. 

Proposition 5.13. If an irreducible billiard-like interval exchange transformation 
f satisfies the modified Keane condition, then f is minimal. 

An example. Consider billiard trajectories within domain T>o with the caustic C ao , 
as described in Section [5. II In addition, suppose the rotation numbers correspond- 
ing to ellipses C\ 1 and C\ 2 are: 

5 1 5 1 



fi = — + — - , r 2 = — 



11 227T 11 220tt 

With given rotation numbers, the Cayley-type conditions from Theorem 13.351 
can be rewritten in a simpler form. Namely, a necessary condition for existence 
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of a trajectory within Vq which becomes closed after n reflections of C\ t and m 
reflections off C\ 2 is: 

Wi + mr 2 G Z. 

In this case, this condition is satisfied for n = 1 and m = 10: 
(5.1) n + 10r 2 = 5. 

1 r 2 

Since — + — — n > 0, the corresponding interval exhange transformation is 
given by: 

(ABCDEFGH 
\DBECHFAG 

/111 21 5 1 111 21 5 I 



V20tt' 2' 22 220tt 11 22tt 20tt 2' 22 220tt 11 22tt 

Proposition 5.14. The transformation (II, A) satisfies the modified Keane condi- 
tion. 

Proof. Suppose that p and p' are two endpoints of the intervals such that p' G' 
{-1,0} and f k (p) = p' for some k > 1. Notice that: 

p = ari + (3r 2 +7^, P' = olr x + l3'r 2 + 7'^ 

for some a, a' G {-1, 0, 1}, /3, /?' G {-1, 0}, 7, 7' G {-2, -1, 0, 1, 2}. 
We have: 

P' = f k (p) =P + km + k 2 r 2 + fc 3 i, 
for some integers k%, k 2 , k% such that k\ + k 2 = k, k\ > 0, k 2 > 0. Thus: 

(5.2) (h+a- a')n + (k 2 + - P')r 2 + (k 3 + 7 - 7') ^ = 0. 
Since r\ and r 2 are irrational, equations (|5.1I) and (|5.2I) must be dependent: 

(5.3) a := fc x + a - a' = -^(k 2 = -^(k 3 + 7 - 7')- 

For each £ <E B U F, either /(£) or / 2 (^) are not in B U F, thus 

(5.4) fc 2 <2fci+2. 

Combining (I5.4j) and (I5.3|) we get 8a < 7. Since k 2 is non-negative, (|5.4p gives that 
a = 0, which leads to fc = fci + k 2 < 3. By direct calculation we check that none of 
the partition interval endpoints is mapped into another one, different from —1 and 
by at most three iterations. □ 



In this example, although the Cayley-type conditon for periodicity is satisfied, 
not only that closed trajectories do not exist, but each of the trajectories densely 
fills the ring between the billiard border and the caustic. 
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Figure 42. Domain T>\ within an ellipse having a "wall" on the y-axis. 



5.7. Unique ergodicity. In this section it will be shown that there is infinitely 
many billiard tables bounded by arcs of confocal conies, such that the corresponding 
flow will not be uniquely ergodic. 

Consider the billiard table T>\ whose boundary consists of ellipse Cp 1 from (|2.1|) 
and segment T = {(0,y) | yjb - #2 < V < y/b-^}, with /3i < fi 2 < b: dVy = 
C/3, U r, see Figure \M 

Fix parameter a : such that (3 2 < «o < b, and take C ao to be a caustic. A 
corresponding billiard trajectory is shown in Figure 23] 




Figure 43. A billiard trajectory in T>\ with an ellipse as caustic. 



Proposition 5.15. Billiard flow within domain T>i and caustic C ao is equivalent 
to the following exchange transformation of the interval [—1, 1): 



(5.5) 



Z + n, f€[-l,-rO 

f + ri-1, f € [-ri,T2 - ri) 

£ + ri, f € [r2 - ri, 1 - n) 

C + n-1, £ e [l-ri,r 2 -ri-\ 

f + ri-2, £e h-n + 1,1), 



1) 
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with ri = p(Pi), r 2 
Proposition \2.1\ 



pifa), and p is the corresponding rotation function, 



Proof. The billiard flow is equivalent to the discrete dynamics of touching points 
of the trajectory with the caustic, with the direction of motion taken into account. 

Introduce metric p on the caustic C ao as in Proposition 12.11 Then, we parame- 
trize C ao by parameters: 

p : C ao -> [0,1), q : C ao -> [-1,0), 

which are natural with respect to p such that p is oriented counterclockwise and q 
clockwise along C ao , and the values p — and q = — 1 correspond to touching points, 
contained in the right half-plane and left-half plane respectively, of tangential lines 
from (0, \/b — Having in mind that reflection on the "wall" T changes the 
orientation of motion, we obtain (|5.5[) . □ 



Map (|5.5p is represented by pair (77, A): 

ABODE 
B E A D C 

A = (1 - ri, r 2 , 1 - r 2 , r 2 , n - r 2 ), 

see also Figure 03] 

A B C 



D 



E 




D 



C 



Figure 44. Transformation corresponding to the billiard within T>\. 

Theorem 5.16. There are billiard tables T>\ and caustics, such that the corre- 
sponding billiard flows are minimal and not uniquely ergodic. 

Proof. The transformation (15.51) corresponds to the Veech example of minimal and 
not uniquely ergodic systems [Veel969j , see also [MT2002] . Namely, choose ao and 
Pi such that r\ = p(/?i) is an irrational number with unbounded partial quotients. 
Then there are irrational numbers r, such that for r 2 = /o(/3 2 ) =r\ — r, the measure 
p on C ao is not ergodic, thus not uniquely ergodic. □ 
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